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Abstract 


If a Space Shuttle flight must be aborted before attaining escape velocity, the 
propellant for .ne payload would be ejected into the stratosphere or the ionosphere 
(which includes the mesosphere and the thermosphere). 

Assuming a linear relationship between the stratospheric loading of NO, and the 
magnitude of the o/one perturbation, we have calculated the change in ozone expected to 
result from the space shuttle ejection of N2O4, based on the ozone change that is 
predicted for the (much greater) NO, input that would accompany large-scale operations 
of SSTs. These calculations show that the effect on ozone is negligibly small. 

The NjO^ may also be released in the ionosphere. Because of the localized and 
transient na:ure of the effects, it is concluded that this will result in no adverse 
environmental impacts. 

We have critically reviewed possible stratospheric fluorine reactions to evaluate the 
magnitude of tluorine-induced ozone destruction relative to the reduction that would be 
caused by addition of an equal amount of chlorine. The predicted effect on stratospheric 
ozone is vanishingly small. 

A similar evaluation was made for an ionosphere injection. No adverse environmental 
impacts are predicted. 
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An Assessment of an F2 or N2O4 Atmospheric 
Injection From an Aborted Space Shuttle Mission 


I. Introduction 

A. ObJ«ctiv«s of the Report 

The purptKe of this report is to evaluate the potential for 
environmental impacts of possible Space Shuttle <SS) payload 
propellants that might be dumped in the stratosphere or iono- 
sphere during a return-to-launch-site (RTLS) aborted mission. 

Reflecting the awareness and concern of scientists and 
engineers for the environment, major technological undertak- 
ings are now preceded by assessments of effects on the envi- 
ronment to allow for planning, rcearch.and design changes to 
mitigate such impacts. Recently, one area of particular con- 
cern has been the ozone-containing stratospheric layer protect- 
ing our planet’s surface from harsh and biologically d> v; 
types of solar radiation (Ref. 1). Some scientists ’oelieve iiie 
ozone layer is threatened by high-flying commercial airliners 
(supersonic transports — SSTs) and chlorofluoromethanes (flu- 
orocarbons), and therefore atmospheric modeling and chemis- 
try have become the subjects of an intensive worldwide re- 
search effort to quantify such effects (Ref. 2). 

While the potential for damage from these sources is still 
being evaluated, the fragile nature of the upper atmosphere has 
become common knowledge (Ref. 3). As such, it was realized 
that large-scale aerospace projects injecting exhaust gases into 


the harth's upper atmosphere presented a potential hazard to 
the natural atmospheric balance (Ref. 4). The intent of this 
investigation has been to estimate the magnitude of this poten- 
tial for one aspect of the SS project. 

Since its developmental stages, various aspects of the SS 
project have been subject to impact assessment. Included 
among these are NO^ generation from SS reentry heat 
(Ref. 5), MCI exhaust gases from SS booster motors (Ref. 6). 
and others (Refs. 4, 7, and 8). 

Tliis report assesses another dimension of the SS project: 
the potential impacts from unburned propellant jettisoned in 
the stratosphere or ionosphere when an SS mission is aborted 
before escape velocity is attained. 

Tlie major portion of this report deals with the chemistry 
and modeling of N2O4 and Fj injections into the stratosphere 
and ionosphere. Both chemical species are treated indepen- 
dently because they represent not a mixture of SS fuel, but 
the principal constituents of two fuels being consideied for 
various payload packages of the SS. While other SS assess- 
ments have considered the case of 60 flights (and injections) 
annually, injections in this report arc modeled as single, inde- 
pendent events because the concern of this report is restricted 
to aborted missions. While it is felt that a series of injections of 
this type could be simulated by a simple “scaling-up," such a 
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irealment is superfluous because any such series of mission 
failures would siMti result in design changes. 

B. Atmosphtric Characttristics 

The Barth’s atmospliere, while consisting of gases, is non- 
uniform in terms of temperature, pressure, and composition. 
Ambient temperature has a range of over 400 K. At the 
surface of the Barth, temperatures average about 2^0 K and 
tlien decrease with altitude due to adiabatic cooling. However, 
at about 17 km, the altitude-temperature '.‘lationship reverses 
atiu the temperature begins to increase with altitude as a result 
of the presence of the oeone layer. In fact, there are three such 
temperature inflections in the atmosphere that occur at ab»ut 
17, 50, and 85 km. Tlie region below 17 km is the tropo- 
sphere. This region is neglected in the remainder of this discus- 
sion because no injections would take place in that atmo- 
sptieric layer. How<r<er. the troposphere is important and 
different from other regions of the atmosphere in that it is wet 
(higli lIjO content). This gives rise to a major atmospheric 
removal process for many impurities by heterogeneous rain- 
out. These temperature intlections. or inversions, result in 
stratification of the atmosphere into Lyers as illustrated in 
Fig. I. Tliis stratification is important in that it results in very 
slow vertical mixing of atmospheric constituents, so that tri-.i 
different layers of the atmosphere exhibit different character- 
istics. Tlie extent of the variation in atmospheric constituents 
with respect to altitude can be seen in Fig. 



AVERAGE TEMPERATURE, X 
Fig. 1. TherniRl profile of the atmoephere 

Not only do the layers differ in lemperalurc and specific 
atmospheric constituents, they also differ in the type of chem- 
istry that predominates. Tlie important chemistry in the 
stratosphere is homogenous, neutral chemistry i rrninated by 
radical species of the HO^. NO^.andO^ families. Perhips the 
most significant constituent of the stratosphere is o/one, O,. 
which shields the Barth’s surface from ultraviolet (UV) light. 


T)ir maximum concentration of atmospheric o/one is observed 
in the stratospliere at about 25 km. 

Above the stratospliere arc the mesosphere and thermo- 
sphere, as seen in Fig. I . This region above 60 km, also collec- 
tively called the ionosphere, is usually further subdivided into 
three regions. The D region is from 60 to 90 km, and ioniza- 
tion there results primarily from the photoioni/ation of NO 
( NO + hr^ -* NO’ ♦ e ). Between 90 and 1 20 km is the B region, 
where photoionization of O, is the most important ion- 
forming step (Oj + hi'-»0* + c). In the F region, above 
1 20 km, O. Oj . and Nj all photoionize, (O ♦ Im -» O’ + e, and 
Nj ♦ In'-* Nj + e)(Ref. 9). 

Tlie approaches used in this assessment reflect the major 
concerns for the different atmospheric layers, which are, in 
turn, a function of the dominant characteristics of the atmo- 
spheric region. The important characteristic of the strato- 
spliere, for the present purpose, is the ozone content: there- 
fore, the principal direction of the analysis for the strato- 
splieric injection is the impact on the ozone concentration. 
Analogously, the assessment of the ionospheric region is prin- 
cipally concerned with the effects of a propellant release on 
the concentrations of thr dominant ions. C'eitain unique fea- 
tures of high-altitude neutral chemistry were alst) considered. 


II. Injection of the Propellant 

A. Infection Parameters 

Tlie opportunities for fuel dumping within the lliglit plan 
for an RTFS abort mode an limited by competing operations 
of higlier priority, as well as factors affecting the safety of the 
crew and vehicle. Such considerations restrict times of dump- 
ing to 153 to 363 s. 442 to 654 s and 744 to 984 s after 
liftoff, as summarized in Table I. These three time spans are 
designated Segments I. II, and III. In the final RTFS llight 
plan, all three llight segments may be retained as alternates, or 
one specific initiation time may be chosen. In this report, 
propellant dumping isevaluaied for all three lliglit segments as 
well as for dump durations of 50 s. 100 s. and 200 s. 

Representative initial plumes that would result from these 
candidate initiation-duration combinations are schenvincally 
depicted in Fig. 3. a graph the SS altitude as a function of 
ground-elapsed time (GBT). The principal effects of the dump- 
ing duration are those of plume volume and initial dilution of 
the propellant, as well as definition of the layers of the 
atmosphere that the plume traverses. Velocity of the vehicle 
has a similar effect on initial concentrations. Trajectory guide- 
lines are summarized in Table 2. Additional detail is available 
elsewhere (Ref. 10). 
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ALTITUDt, km 



Fig. 2. Concentration* of chemical apeciet In the atmoaphere according to altitud* 


Table 1. RTLS flight plan c^gmenta aatiafying criteria 
for propellant dumping 


1 lieht 
SeeiiK'iit 

RTLS 

events 

Gi r, 

s 

.\ltitude. 
km tt't) 

1 

Time of abort 
to pitch around 

153.9 to 

363.9 

63.3 to 135.3 
1303.976 10 411.188) 

II 

Masinium ran):e 
distance to 
Ml CO 

443 to 
654 

91 4 to 70.0 
(300.000 to 339.756) 

III 

X aeeleralion 
control to 
Mach No. 3 

744 to 
9H4 

43.9 to 38.5 
(140.890 to 93.458) 


The two partieulur SS payload propellanls being evaluated 
eontain 2722 kg (h.OOOlb) of I ,, or .'402 kg (7.5tX)lb) of 
S,0^. I or either propellant, this is a worst ease analysis with 
regard to quanttty iti that lesser, but not greater, antounts of 
these substances may eventually be carried by the SS payload. 

In terms of relevant aspects of the injection altitudes, all of 
the itiitiation-duration dutnpitig optiotis may be generali/ed as 
two cases. Dumps in Segments I and II would be restricted 
entirely to the ionospheric layer of the atmospheie. which 
includes [xirtions of the mesosphere and the thermosphere. 
,\ny injections in Segtiient III would be restricted to the 
stratospheric layer. Hence, the two injection cases evaluated 
wete an ionospheric injection and a stratospheric injection. 
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Fig. 3. Tr a ttc t ory of Iho Spaco SEmiMo during an RTL8 abort 
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Those chemical species conceivably decreased by the 
injected propellant include stratospheric o^one and iono- 
spheric charged particles. The magnitude ot any perturbation 
to the ambient atmosphere would be dependent upon the 
concentration of the pollutant, which is a function of the 
absolute amount of the injected material and the rate of 
dilution due to transport processes. After a period of weeks to 
months. ste;.Jv state (equilibrium) conditions would be 
reached and 'i.e injectant will have hemispherically dispersed. 
During this same time period, very little vertical transport will 
have occurred. Therefore, the steady-state concentration of 
the propellant merely reflects the heigitt of the altitude regime 
into which the injection took place. 

On a shorter time scale, before the plume hemispherically 
disperses, pollutant concentrations will be substantially higher, 
as will be the accompanying transient and localized perturba- 
tions. An evaluation of atmospheric diffusion parameters has 
been made in an attempt to estimate the time scale of the 
plume dispersal and the intermediate propellant concen- 
trations. 

B. Plume Dynamics 

The main effects of fluid motions on atmospheric injectants 
arc zonal (E-W) and meridional (N-S) circulatory advection. 


and turbulent eddy diffusion. Tlte primary effect of the circu- 
latory advection is simple horizontal displacement, whereas 
turbulent eddy diffusion causes dispersion of the injectant. In 
the stratosphere and mesospitere, the lurrizontal velocities (~1 
to 100 ms~' ) are one to six orders of magnitude greater than 
the vertical velocities (~ 10 "* to 10 "' ms*'), and arc not 
constant at different altitudes with respect to velocity or 
direction (Ref. 1 1 ). The wind parameters abruptly change as is 
seen in Fig. 4 (Ref. 12). The effect of these abrupt changes is 
one of vertical stratification. Dispersion of the injectant is 
produced by both small and large scale eddies (in the strato- 
sphere and mesosphere), and can be parameterized in terms of 
an overall eddy transport coefficient (A) by assuming that the 
eddy tlux (f'„) is proportional to the gradient of the mean 
concentration of the pollutant being diffused (V/i^). such that 

Fa = -KV,la 

Tlicorctical values of the horizontal (A^,^, ) and vertical (A..) 
eddy diffusion coefficients are available in the literature, based 
on the observed distributions of trace gases such as O 3 . IKO. 
N,0. CH 4 . and radioactive debris. Figure 5 shows typical 
values of Kyy and A'., generated in a two-dimensional model 
(Ref. 1 1 ). It can be seen from these values of A (~3 X 10'* 
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TaiiM 2. R tlM aboft; avant aaquancia dala 


Pvem 

GET, « 

Rebtive 
velcKity 
y, km/t(ft/i) 

Altitude 
h, km (ft) 

Time of aburr 

1S2.9 

1.38(4515) 

62.17(203.967) 

Tu.n around 

362.9 

2.55 (83501 

125.33(411,188) 

i-g't accekialion 

622.9 

1.43(4706) 

69.67(228.588) 

Pilch down 

649.9 

2.07(6781) 

69.40(227.676) 

Main eniune cutoff 

6S4.4 

2.19(7193) 

69.81 (229.022) 

(MfCX)) 

External tank lep- 

661.4 

2.26(7405) 

70.03(229.756) 

aration (ET SEP) 

ET SE P 

661 

2.26(7405) 

70.10(229,998) 

Pitch up 

671 

2.26(7410) 

69.61 (228,383) 

Angle of attack, a/f 

691 

2.27(7449) 

65.85(216.047) 

Dynamk' prruure. 

712 

2.27(7463) 

58.43(191.669) 

</ • 20 

Mach No. » 7 

716 

2.27(7459) 

56.66(185.891) 

K acceleration control 

744 

2.16(7092) 

42.94(140.890) 

Mach No. « 6 

771 

1.85 (6058) 

34.95(114.680) 

.Maximum q 

774 

1.81 (5946) 

34.65(113.689) 

Pull out (•) ■■ 0) 

781 

1.73(5681) 

34.37(112.771) 

Mach No. » 5 

799 

1.55 (5081) 

35.79(117.411) 

Mach No. • 4 

883 

1.23(4032) 

34.03(111.634) 

Mach No. = 3 

948 

0.90(2950) 

28.48 (93.453) 

Terminal augmented 

1054 

0.46(1499) 

21.24(69,689) 

energy manage* 
men! (TAtM) 

Autolanding 

1375 

0.17(544) 

2.67(8,754) 

Landing 

1441 

0.09(295) 




cm*-» ' ) anti A',, (~l to 20 X 10^ ) tlut huri/ontal 

niixtng » significantly more rapid than vertical mixing. How- 
ever. the actual diffusion coefficient values used in atmo- 
spheric models vary considerably, as is evident from Pig. b, a 
graph of A'j coefficients used by different groups in their 
one-dimensional models. Figure 6 also illustrates the variation 
in A'j with altitude. 

The vertical diffusion cirefficienl. A,, can be used to calcu- 
late the residence time, which is an indication of how lor^! 
pollutant will remain m the atmosphere. The residence \-,ni 
defined as the time it takes for the amount of the pollutant 
that remains in the alm.osphere to decrease to I/e of the 
original amount (the e-folding tiiiie). These calculations nor- 
mally assume that the lifetime in the troposphere is sliort (a 
few days) and controlled by heterogeneous rain-out priKesses. 
Figure 7 is a graph of the residence times calculated from the 
profiles shown for A^ in Fig. 6. It can be seen that tlie 
predicted residence times in the Chang and McElroy, et al., 
models (Kef. 1.^) differ by about 8 years for an injection at 
45 km. This illustrates the importance of the A^ profile 
chosen. If the McFlroy, et al., profile is used, there is essen- 
tially no difference in residence time between a 22-km injec- 
tion and a 45-km injection. However. A. values used in the 
Chang model result in a difference in residence time of over 10 
years for those same injections. Relative to the McFlroy, et al., 
model, the Chung model would predict a longer residence time 
for mesospheric injections. Recently, however, the consensus 
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Rg. 5. Th« variation in value* tor transport coelficlenta with 
aKitude and longitude for March through May: (a) horizontal 
(K^) distribution (10^ cm^ (b) vertical (K„) distribution 
(10^ cm2 ,-1) 

among the modeling eommunity is lliat ihe best A', profile is 
similar to that used by Mcblroy. el 

I'igiires 8 and are graphs of the A', profile and resultant 
residence times used by Yung (Ref. '4( in his model of the 
atmosphere. We used ^'ting’s model to obtain the concentra- 
tion profiles in f'igs. 10 and 11. These profiles illustrate the 
deciease in concentrations and mean altitude of the pollutant 
following a 40-km and a l(X)-km injection. Based on this 
model, a pollutant injected at 100 km docs not reside signifi- 
cantly longer in the atmosphere than for a 40-km injection. 
Since the eddy diffusiim coefficients used in the Yung model 
are similar to those of Mctlroy. el al.. bi th predict that the 
residence time is approximately constant above 25 km. It is 
probable that Ihe model (Ref. 141 used to obtain the concen- 



Fig. 6. profilM UMd by diffarvnt atmospharic 
ipodaling groups 



RESIDENCe TIME t, yr 

Fig. 7. Rssidertes time profilet calculated using Kj coefficients 
from Fig. 6 


liaiion profiles in Figs. 10 and 1 ' indicates a more rapid 
vertical dispersion than w'ould be predicted by other models. 
Wliile these differences are of academic interest, the resulting 
impacts calculated in this evaluation would not appreciably 
change if an alternate A', profile were used. Figure 12 illus- 
trates h.ow the mixing ratio of a pollutant injected at 100 km 
decreases with time, due to ver'ical dispersion, and becomes 
more uniform with altitude. 

Above I(X) km. Ihe total particle density decreases and 
molecular diffusion becomes the dominant mechanism for 
vertical mixing as illustrated in Fig. 13 (Ref. 15). Tlie rate of 
dispersion in this region of the atmosphere can be described by 
the equation lor spherically symmetric diffusion (Ref. 1(>); 

OF POOR 









Fig. 9. Residence tinie profiie calcuiated using coefficients 
from Fig. 8 



where; 

n = particle density 
r = radius of the cloud 
t = time 


Tlie solution tan be written: 

»(r,f) = - — - exp (-r^/4Dr) 

{4nDl)^'^ 

where 5^, is the total number of molecules injected and D is 
the molecular diffusi(’n coefficient. This mathematical treat- 
ment of plume dispersion was not used in the present report 
because no propellant, or only a very small portion of the 
propellant, would be injected above 100 km. 

C. Plume Disperse! in the Stratosphere 

The first situation to be consiilercd is that where the SS 
mission is aborted in the stratosphere. The propellant can be 
released in any 50-s. lOU-s, or 200-s time period, while the SS 
is in the altitude range of 43 to 28.4 km. Tlie trajectory of the 
SS during the injection mode and the duratior of the dumping 
procedure determines the altitude regime of the injection. As 
the rate at which the plume disperses is predominantly a 
function of rapid horizontal mixing processes, the rate of 
plume expansion is enhanced by deposition over a wider 
altitude range. 

Four different stratospheric propellant releases are consid- 
ered; 50 s. 100 s, and 200 s releases initiated at 744 s (itT (the 
earliest opportunity in Segment 111), and a 50-s release initi- 
ated at 771 s GtT. The releases initiated at 744 s correspond 
to an altitude ‘d’ 43 km. and the 771-s release occurs at an 
altitude of 37.6 km. Tlie release initiated at 771 s corresponds 
to the situation where the propellant is deposited in the 
narrowest possible regime, when the SS is on a near-horizontal 
night trajectory. 

Table 3 lists values for the critical parameters involved in 
the four cases; altitude regime, distance travelled by the SS 
during 'lie ejection period, ejection duration, anu initiation 
time (GFD. In each of the cases, (a) througli (d). the amounts 
of Fj and NjO^ are 4.2 X 10“ and 2.2 X 10“ molec iles. 

The N 2 O 4 rapidly undergoes thermal decomposition lesulting 
in the productio.T of 4.4 X 10** molecules of NO,. 

Aithougli the diameter of the propellant exit port is postu- 
lated at only 2.5 cm. the plume diameter rapidly (less than one 
second) expands by a factor of 50 to 200 due to the pressiiu- 
differential between the ejected p’unie and the ambient atmo- 
sphere In case (a), the gas pressure in a 2.5-cm diameter plume 
would be ~8500 times the total pressure of the ambient 
atmosphere. In cases (a) and (d) the initial plume diameter 
would be ~ i to ' i. whereas in cases (b) and (c) the plume 
diameters would 2 to 2-1/2 and ~1 m respectively. Tlius 
in cases (b) and (c) the plume does not expand as much, due S 
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Fig. 10. ProfllM of propoilant concontrations foilowing an Iniaction at 40 km 



Fig. 11. Profiiet of propeiiant concantrationa foiiowing an injaction at 100 tun 

8 


obigwalpa^ 

OF POOR QUAUr* 







CeOMEnttC ALTITUOe, bn 



Fig. 12. Mixing ratio of tha ln|actant at diffarant altitudaa aa a function of tima following an inioction at 100 km 



Fig. 13. Dominant machanisms inducing pluma dlaparsal 
according to aHKuda 


to a lower pressure differential between the plume and the 
ambient atmosphere. 

As soon as the gas is ejected into the atmosphere, it under- 
goes shearing by horizontal winds and dispersion due to turbu- 
lent diffusion. Figure 14 illustrates the horizontal and vertical 
mean “cloud” widths, and turbulent diffusion coefficients. 


Tabla 3. Critical paramatara for aampia 
atratoapiiaric Injactiona 


Case 

Altitude, km 

Duration, s 

Disianee. km 

Initiation times 

(a) 

43.0 to 35.7 

50 

92.9 

744 

(b) 

43.0 to 34.4 

too 

175.8 

744 

(c) 

43.0 to 28.4 

200 

.’oe.! 

744 

(d) 

37.6 to 38.5 

50 

ec 5 

771 


Kyy, Ky,, and A'jj. as a function of travel time (plume age) 
for a hypothetical point source injection (Ref. 17). It cue. be 
seen from this figure how the magnitude of the horizon- 
tal eddy diffusion coefficient, varies with time (or plume 
diameter). At short times the dispersion is isotropic and ap- 
pro dies the molecular limit fur the injection altitude, and 
ihen Kyy increases as o}, ' * up to o^, ^ 100 km (Of, = horizon- 
tal dispersion cloud width), at which point it reaches the 
global limit. Tlie global limit shown for A' in Fig. 14 corre- 
sponds to a value of ~3 X 10*® cm^-s"', which is approxi- 
mately an order of magnitude greater than the values calcu- 
lated by Louis (Ref. II), but comparable to those calculated 
by Crutzen (Ref. 18). in their two-dimensional models. The 
values of Kyy calculated by Louis and Crutzen were typically 
3X 10® cm2-s'*,and > 10'® cm^-s”* respectively, and repre- 
sented meridional dispersion (N-S). Zonal dispersion is prob- 
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ably comparable in magnilu te. Although Fig. 14 has been used 
to estijnaie the rate of '.lori/.ontal dispersion, it has been 
realized that the data base used to develop this figure is limited 
at best. If Fig. 14 is in error at short times (r< 10* s). it 
would probably be on the side of underestimating the rate of 
plume dispersal due to an incorrect time evolution for 
However, at longer times. Fig. 14 may overestimate the rate of 
dispersal if the global asymptotic value for Ky^, is too high. 
For the region where A',,,, or A'., have reached their asymp- 
totic values, the rate of plume dispersal can be calculated using 
the following simple expression: 

o' (mean plume width = 4 A / 

Figure 14 has been used in this assessment to calculate the 
time history of the plume volume. In this figure and our 
subsequent calculations. A'^,^. represents the parameteri/ed 
eddy diffusion ciK'fficient which describes both small and large 
scale eddies. It ran he seen J'rotn I 'ig. 14 that the plume size at 
lnnf! times is independent of any assumptions made about the 
initial shape and size of the plume. Whether the initial cloud 
width is 1, 10. or 100 m, it will still expand to '-10- km after 
1 day and ~10* km after 10 days. Consequently, the esti- 
mated plume width after 1 day, 1 month, and 1 year were 


taken di: ctly from Fig. 14. \t short times (< I day) the 
calculated pluiiic volumes may be somewhat dependent upon 
the validity of various approximations. As stated earlier, the 
initial diameter of the gas plume (within a few seconds after 
release) is dependent upon the ratio of the “plume pressure” 
and the ambient pressure of the atmosphere. The plume will 
rapidly expand to eliminate the pressure gradient. In addition, 
as the SS takes either 50. 100, or 200 s to release its propel- 
lant. the gas released at the beginning of the period has more 
time to disperse. Consequently, the shape of the plume at the 
end of the period can best be described as the tilted fmstrum 
of an inverted cone (see Fig. 1 5a). 

Althougli the SS travels ^^3 km during this period, it 
descends only ~7.3 km in altitude, largely cimfining the prr>- 
pellant to this altitude regime for several months due to the 
slow rate of vertical mixing. Consequently the plume is pre- 
dicted to e.x/tand rapidly in volume by horizontal dispersion, 
hut only very slowly by vertical mixing. Tltis merely reflects 
the large difference in magnitudes between the rates of lurri- 
zontal and vertical mixing. Within a few days, the shape of the 
plume changes from the tilted frustrurn of a cone (shown in 
Fig. 1 5a I to that of a cylinder illustrated in Fig. 15b. For times 
between the end of the dumping period and the time it takes 
for the plume to reach the shape of a cylinder, the plume 



Fig. 14. Dispersion distances and coefficients of turbulent diffusion as a function of time following an injection at 100 km 
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Fig. 15. Qconwtric shape of the propsllant piums: (a) Initial 
shape of plume, (b) plume shape after a faw days 

shape and rate uf expansion are less easily calculated. How- 
ever, it is quite adequate and intuitively correct to assume that 
the shape of the plume can be described, for the purpose of 
calculating its rate of expansion, by a cylinder ~7.3 km in 
heiglit with an initial volume equivalent to that of the initial 
tilted inverted cone. 


Tliese calculations of plume volume are probably good to 
within only an order of magnitude because of uncertainty in 
Kyy. Hie plume dispersal rates will vary greatly depending 
upon the meteorological conditions prevailing at the time. Kyy 
is predicted to exhibit significant variability with both latitude 
and season (Kef. II). Table 4 shows the calculated plume 
volumes and resultant FX and NO^ concentrations and mixing 
ratios as a function of time, assuming no vertical mixing. This 
approach is obviously invalid at times greatei than I month 
and leads to a small underestimate of the true plume volume. 

0. Plume Oisperul in the Mesosphere 

Tlie propellant can be released during any 50-. I00-, or 
2(X)-s lime period while the Space Sliutile ascends from 
62.2 km to I25.‘)km (I.S2.‘)— .362.6s Cil-T), or descends 
from 91 .4 km to 70.0 km (442 — 654 s GFT). Any 50-s dump 
in night Segment I will initially be distributed within an atmo- 
spheric layer 15 km in heiglit, and over a hori/onlal distance 
of ~79 km. Tlie initial volume of the plume (at the end of the 
dump peiiod) is governed by the almost instantaneous expan- 
sion that iK'curs due to the pressure differential between the 
ejected propellant and the ambient atmosphere (the initial 
expansion is greatest at higher altitudes due to the lower 
atmospheric pressure). Let us consider two extiei.>e cases: 
(l)a 50-s dump between 62 and 77 km, and (2) a 5C-s dump 
between III and 126km. It is easily shown that the initial 
diameter of the propellant plume will be '-42 m in case ( I ) 
and 2.54 km in case (2). From Fig. 14. it can be seen that 
there will be no significant exp ;nsion of the plume within the 
dump period due to eddy processes. (This is in contrast to the 
stratospheric case.) The initial plume volumes are .3 X 10' -’ 


TabI* 4. Tima history of ptiime volumes and concentrations 
of FX and NO, for a stratospheric injection 


Time 

Plume volume, 
cm^ 

Cloud width . 
km 

|I■■.X|“, moleeules/cm^ 

|NO^ 1 , moleculcs/cm^ 

(1 X)* 
ratio 

nuxing 

ratio 

0 

3.2 X 10*^ 

2.4 X l()-3 

1.34 X lO'* 

1.38 X 11)** 

0.17 

0.17 

1 h 

1.0 X 10** 

1.4 X 10“ 

5.0 X 10*3 

4.4 

X 10*3 

6.3 X 1()-* 

5.5 X 10** 

5.6 h 

1.0 X 10*® 

1.4 X 10* 

7.1 X 10*° 

4.4 

X 1()*° 

8.9 X K)-^ 

5.5 X 10*'' 

I day 

5.1 X lO'** 

9.6 X 10* 

1.66 X 10** 

8.6 

X 10* 

2.1 X !()-* 

l.l X It)-* 

10 days 

8.7 X 10*' 

1.26 X 103 

9.8 X 10* 

5.1 

X 10* 

1.2 X l()-*° 

6.4 X itr* * 

1 mo 

7.4 X 10^3 

3.7 X 103 

1.16 X 10* 

5.9 

X 10* 

1.4 X 1()-** 

7.4 X 10** 3 

1 year 

1.8 X 103“* 

Hemisperic 

4.7 X lO"* 

2.3 

X in'* 

5.9 X !()-** 

3.1 X 10** •* 

5 years 

3.6 X 103“* 

Cilobal 

2.3 X 10“’ 

1.25 X 10“* 

2.9 X 10*** 

1.6 X 10**-* 


“I X = III I + I O FOj + HOI + NOI 

*’lt is assumed that the cloud center remains al 40 km. there is no vertical dispersion, and the propellant is uniformly mixed within the cloud. Both 
assumptions are incorrect, but arc used for simplicity (discussed with reference to Fig. 11). |M| = 8.0 X 10** molecule-cm"-* (total gas content 
al 40 km). At times less than one day. the values shown for I X have taken into account the I'j photo dissociation rale (each molecule of I'j pro- 
duces 2 l-'X molecules). 
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cnv* und 2 X 10' *• cm-' rcspcLtively. Tlie plume width at times 
lunger than several hours is nut dependent upmi the initial si/.e 
of the plume (Hg. 14). Consequently, the time history of the 
cloud width will be the same for any dumping altitude or 
dump duration. The volume is dependent on the dump dura- 
tion as this governs the initial heiglit of the plume. Table S 
shows the predicted time history o*' the plume width, plume 
volume. |h'X|,and |NO^|. 

TaM* S. Ptunw vohimat and conca'itratkma ol FX and NO. 
for an in)acttcn In tha lonoapfiara 


Plume width. Plume volume, 
T ime . 1 

km em' 


1 h 

2 

4.7 

X 

lo'*’ 

1.8 

X 

10'2 

9.4 

X 

10" 

1 day 

96 

I.l 

X 

10^' 

7.8 

X 

10’ 

4.0 

X 

lO’ 

10 days 

1.26 XlO-* 

1.9 

X 

10“ 

4.6 

X 

10‘ 

2.3 

X 

6 

10 

1 mo 

3.7 XIO^ 

1.6 

X 

10^^ 

5.4 

X 

10* 

2.8 

X 

10* 

1 year 

Hemispheric 

3.9 

X 

10“ 

2.2 

X 

10^ 

1.1 

X 

10« 


From these values of plume volume and |F.X| (it has been 
assumed 2 FX molecules arc produced per F^ injected), the 
mixing ratios (r/r) of FX at different injection altitudes were 
calculated (Table 6). When the propellant is released over a 
100- or 200-s lime period, these values must be reduced by 
factors of 2 and 4 respectively. In addition, these values 
neglect vertical dispersion, which is quite significant at times 
greater than I month at these high altitudes. Figures 1 1 and 12 
show how vertical plume dispersion is quite rapid for higli 
altitude injections. Therefore these values tend to overesti- 
mate the mixing ratios of FX and N()^ (almost identical to 
FX - see Table 5). 


Table 6. The mixing ratios (v/v) of FX at different altitudes 
following an ionospheric injection 


Time 

70 km 

90 km 

120 km 

I day 

3.9 X 10’* 

1.0 X 10"* 

1.4 X lO-* 

10 days 

2.3 X 10"’ 

6.0 X lO"" 

8.5 X 10*® 

1 mo 

2.7 X 10‘'“ 

7.0 X 10*’' 

1.0 X 10"* 

1 year 

U X ur" 

3.0 X I0*'° 

4.1 X ur* 


III. Atmospheric Modeling 

A. Background 

Basically, a stratospheric model is a set of equations 
representing all constant and dynamic aspects of the strato- 
sphere that have a direct or indirect effect on the concentra- 


tions of the relevant chemical species. Representative mathe- 
matical expressions are used to generate a computer program. 
Tliese models have become increasingly complex, to the point 
that a run may take hours or days of computer time. 
Typically, a model takes years to develop and will require the 
concentrations of thirty chemical species while utilizing rate 
constants, k and J values (for photiK'hemical reactions), for 
100 or more reactions. Additionally, the models must take 
into account natural variations in oackground levels, diurnal 
variations of sunliglit, the horizontal and vertical movements 
in the stratosphere, and the altitude and location of the 
pollutant injection. Finally, many of the concentrations, rate 
constants, and other input data are unavailable and must be 
estimated. Thus, it is easy to see that two groups of modelers 
could make different estimates and predict different out- 
comes. In fact, most recently, not only the magnitude, but 
also the sign of the impact of the SSTs is in dispute. S.ime 
groups presently claim '.!.it SSTs flying at 17 km would 
increase rather than decrease the level of stratospheric zone 
(Ref. 19) 

Models exhibit varying levels of sophistication or 
"parameterization" according to the level of data-averaging 
that is involved. The least complex models are 1-D( dimension) 
models, or vertical profiles. Tliese models emphasize complete- 
ness with regard to the chemistry of the stratosphere and are 
very comprehensive with regard to the reactions in the model. 
Vertical transport is also emphasized. However, horizontal 
movement, on a micro or macro level, is ignored and 1-D 
models must assume a globally-averaged input. Additionally, 
vertical transport as well as the ambient concentrations must 
be averaged to account fur geographical variations. Ihus. the 
1-1) imxlel cannot predict the extent of local or “corridor" 
perturbations (disruplions) in the stratospheric chemistry 
because the output is only a vertical profile. However, the 1-D 
models use a fine vertical grid size, thereby allowing for more 
detail in the vertical direction. They also allow for more 
completeness with regard to photochemistry. 

The 2-D models represent a cross-section of the strato- 
sphere througli a polar axis. While 2-1) models can incorporate 
more detail for K^y (the zonally averaged meridiona! eddy 
diffusion coefficient representing north-south horizontal trans- 
port). often fewer reactions arc included and the vertical grid 
size may be larger. Tliis addition of Kyy is a substantial 
addition in that interhcmispheric transport is relatively slow 
and the altitude of the ozone layer exhibits variation according 
to latitude. 

Finally, there are 3-D models. These models make maxi- 
mum use of horizontal diffusion in both the polar and 
equatorial (A\) directions. However, still more detail in the 
chemistry and subdivision of the stratospheric layer must be 
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sacrificed. Furthermore, the demand for computer time to run 
the model substantially increases. 

B. Assessing an NOi ln)ection 

The model for NO, injections has undergone intensive 
development under the Climatic Impact Assessment Program 
(CIAPI. Department of Transportation. Under CIAP, several 
nio'icls have been developed to predict the impacts of present 
subsonic and future levels of SST fliglit on the stratosphere. 
Thus, the NO, models are well establislied and. while some 
areas of controversy still remain over effects of aircraft 
injections at lower stratospheric levels, the models are in 
substantial agreement regarding the overall role of NO,. (The 
agreement between I-. 2-. and .TD models is not unexpected in 
that, althougli lower-order models use averaged data for more 
of the parameters, the data base itself is the same for all of the 
models. I 

In this assessment we have taken advantage of the advanced 
state of NO, models. In so doing, we have used the results 
from existing models to predict the ozone perturbation, at 
steady state (equilibrium) conditions, expected to result from 
an NO, injection from the SS payload. This analysis is 
presented m Section IV. 

C. Assessing an FOx Injection 

In the course of assessing the impacts of the MCI injected 
into the stratosphere by the SS main engines (not to be 
confused with the present case of a SS payload propellant), 
models were developed for stratospheric chlorine chemistry. 
From that HCI assessment, the catalytic nature of chlorine 
destruction of ozone was realized, even though the HCI 
injected by the main engines appeared not to pose a serious 
threat to the ozone layer (worst case predictions are '•0.5% 
decrease in ozone) (Ref. 6). Rowland and Molina (Ref. 20) did 
a study indicating that chlorofluoromethanes (CFMs) can 
pliotodissiK'iate in the upper stratosphere generating Cl. but 
not in the troposphere. Tliis causes long atmospheric lifetimes 
for the CFMs (to date no tropospheric chemical sinks have 
been identified), allowing long-term stratospheric accumula- 
tion of Cl. Due to the tremendous quantities of CFMs used as 
aerosol propellants. CFMs were realized as a potentially very 
substantial threat to the ozone layer. Hence, much like NO,, 
models have been developed for stratospheric CIO, chemistry 
and for stratospheric injections of Cl. Tlicse CIO, models arc 
also well established. WliiJe controversies remain over details of 
the models, the predictions of the various CIO, models are 
substantially the same (Ref. 21 ). 

Althougli chlorine and tluorine are different substances, 
they are of the same family, i.e.. halogens, and exhibit very 
similar chemistries. In assessing the impacts of the lluorine 


injection, we have taken advantage of the advanced state of 
CIO, models as well as the analogous qualities of Cl and F 
chemistries. 

It is not realistic to merely substitute the value for the 
tluorine injection into a chlorine model. Rather, stratospheric 
tluorine chemistry has been very carefully analyz'd to 
generate analytic expressions for the catalytic eft'iciency of 
tluorine comparable to those of Cl. Using experimental and 
estimated rate constant data, the ratio of the values for the F 
and Cl catalytic etTiciencics has been calculated. This factor is 
then used to convert tiie ozone perturbation predicted for a 
chlorine injection in a CIO, mode, to that predicted for a FO, 
injection of the same magnitude. Thus, the full benefit o'' a 
sophisticated stratospheric model can be realized in detem.in- 
ing the ozone perturbation resulting from the fluorine 
injection. 

Details reg. ding analytic expressions, rate constant data, 
and the tluorine model, as well as ihe actual quantities 
calculated and overall assessment, are presented in Section V. 

IV. N2O4 in the Atmosphere 

A. Background 

An NO, model has been developed for higli flying aircraft 
and is fairly well established. Figure I (1 is a schematic diagram 
of the chemistiv use I in the NO, model. Tlie principal 
chemical reactions incorporated into the NO, model are 
summarized in Table 7. 

Althougli the pollutant is injected in the fonn of N, 04 . 
reaction (18). 

N,(>^ -I- M -»2NOj -t- M (18) 

is very rapid under stratospheric conditions and the entire 
amount will dissociate within I to 2 seconds to form NO,. 
Tlierefore. for purposes of analysis, the dumping can be 
considered an NO, injection, the form similar to aircraft 
exhaust. Furthermore, since NO and NO, cause the catalytic 
destruction of ozone througli reactions (II) and ( 1 2), they are 
both considered to be in the “active" partition and -ire 
collectively termed NO,. The basic mechanism for production 
of ozone in the stratosphere is as follows (Chapman 
mechanism ): 

0,-t^ In' -* O + O ( I ) 

O -I- Oj M Oj M. twice (2) 

Net: 30 j hi^ -* O^ -i- 
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t able 7. Stratospheric NO, model 


Reaction 

i'.quatiun 

No. 

Oj + hi’ 

-0 + 0 

(1) 

O + Oj + M 

- O3 + M 

(2) 

O3+O 

- 0, + 0, 

(31 

Oj + ill’ 

- Oj +0 

14 ) 

0 ('D) + lIjO 

- 110 + HO 

( 5 ) 

HO + Oj 

- HOO + Oj 

(6) 

IKK) + Oj 

- HO + Oj + Oj 

( 7 ) 

HO + O 

- H + Oj 

(8) 

1100 + 0 

- HO + Oj 

( 9 ) 

H+O3 

- HO + Oj 

(10) 

NO + O3 

- NOj + Oj 

ill) 

NO, + 0 

- NO + Oj 

(12) 

NO, + ill’ 

- NO + O 

( 13 ) 

HO + NO, + M 

- HNOj + M 

( 14 ) 

UNO, + hi’ 

- HO + NOj 

(IS) 

HOO + NO 

- HO + NOj 

( 16 ) 

NO, + O3 

- NO3+OJ 

( 17 ) 


Tills reaction set produces o/one in the stratosphere at a total 
rate of about S.O X 10-*' niolecules/s ( Ref. 22). 

In a pure oxygen system, ozone is destroyed by reaction 
(3). 


O 3 + 0 - Oj + Oj 

(3) 

Altluiugl) ozone can be rapidly photolyzed (4), this does 
not result in the loss of -vone as reaction (4) is usually 
followed by ( 2 ) 

0 

c 

t 

(4) 

0 + Oj + M -+ O 3 + M 

( 2 ) 

Net; null reaction 


lienee, the rate of ozone destruction is a function of the 
relative rates of (3) and (2) with respect to |0| . Tliis net cycle 
(4 followed by 2) is important in controlling the temperature 
stability of the stratosphere due to the exothermicity of 
reaction ( 2 ). 

The 110^ family provides additional reaction-pairs that 
cause the catalytic destruction of ozone; three examples are 
shown ( 6 )- (7). ( 6 ) (^>). and (10) ( 8 ). 

IK) + O 3 - HOO + Oj 

( 6 ) 

IIOO + Oj-HO + Oj +Oj 

(7) 

Net; 2 O 3 -» 3 Oj 


HO + O 3 -* IKK) + Oj 

( 6 ) 

HOO + 0-HO + Oj 

(4) 

Net: O 3 +0-Oj +0, 


11 + 0, -+H0 + 0, 

110 ) 

HO + 0-11 + 0, 

( 8 ) 

Net; 0 ^ +0 -+0j + O, 
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However, these three cycles must compete ‘or H()^ 
zero balance cycle. 

with a net 

HO + Oj-lKKI + O, 

(6) 

HOO + NO -* HO + NO, 

(16) 

NO, + In' - NO + O 

(13) 

O + Oj + M -« O3 + M 

(2) 

Net: null reaction 


A comparable catalytic destruction cycle and competing 
neutral cycle are demonstrable for NO^ chemistry. 

NO + O3 - NOj + O3 

(II) 

NOj +0-N0 + 0j 

(12) 

Net: Oj +0-03+03 


NO + O3 - NO3 + O3 

(11) 

NO3 + In' -<• NO + 0 

(13) 

0 + O3 + M - 0, + M 

(2) 


Net: null reaction 


Thus for these NO^ eyeles the relative rates of reaetions ( 12 ) 
and ( 13 ) determine how imieh ozone will he destroyed. 


B. NOx Modeling 

As is the case with any body of knowledge, NO^ models 
have evolved thrt)ugli a series of changes to reach their present 
form, and continue to evolve as values for rate constants and 
transport coefficients are improved, and as new reactions are 
added to the models. Table 8 and F'ig. 17 summari/.e the 
evolution of the predicted change in stratospheric o/one In 
response to new or improved rate constant data and addition 
of new reactions to the model. It is evident that the most 
recent change, the new rate constant for reaction (161, 

llOO + NO-‘HO + NO, (16) 


TabI* S. Summary of changaa In pradictad oiona lavala for 
tlw atraloapfiaric NO, modal 


Oianiic made 

Computed 
ozone reduction, % 

1 7 km 20 km 

Start 11974) 

4.8 

~I 1 

New- Chanii iRcf. 23) eddy uiffuiiun 

5.4 

12 

(N 3 O ♦ In') changed from 
Bates and Hayes to I Kef. 24) 
Johnston .^nd Selwyn iRef. 25) 

5.3 

12 

Add "smog" reactions: OH CO, 

0 + CH 4 . HjOj ♦ 0 , HOj ♦ hv. OH ♦ CH 4 

4.8 

11 

0H + H02(2x I0*'“-2x lO"'*) 

2.1 

6.4 

OH + HNO,(I.3x 10 '^-8.9x lO*''*) 

2.0 

6.0 

NOj + hi' (branching. NO + Oj -• 
2/3(NOj +0). l/3(NO + Oj)) 

1.6 

5.2 

0 ( * D) reactions with N 2 O, Nj, 

Oj, CII 4 , II 2 O changed from 
Hampson and Ciarvin ( Ref. 26) recommenda- 
tions to Streit ct al.. (Ref. 27) 
recommendations 

1.5 

4.9 

OH + NO 2 + M (Tsang (Ref. 28) — 
Anastasi (Ref. 29) et al.) 

1.5 

4.8 

IIO 2 + IIO 2 (3 X 10*' ‘ exp(-500/T) 
- 1.7 X 10-11 e\p(-500/T)) 

1.2 

4.3 

H 02 + 0 ( 8 X 10‘" e\p(-500/T) 
-3x !()-")) 

1.2 

4.2 

IIO 2 + NO (Howard and I vcnson. Ref. 30) 

+ 1.9“ 

+.4“ 

“Ozone increase. 



Fig. 17. Evolution of the predicted change in stratospheric 
ozone from stratospheric SST ’light. 
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has I'haniscil not only the magnitude, but also the direction ol 
the o^iAie perturbation fur 17- and 20-krn injections; i.e., the 
models now predict an increase in the integrated o/.une 
column rather than a decrease (Kef. 31 1. 

Another recent development in modeling theory was an 
uncertainty analysis whereby an attempt wis made to quanti- 
tate variability in the models that is due •.<: mcertainties in the 
available rate-consiant data. In such an .'•nalysis, Stularski, et 
al.. (Ref. 32) assigned a normal (Ciaussian) distribution to the 
levels of uncertainty for reaction rate constants. Performing 
several hundred model runs while individually varying rates for 
each of the 48 reactions used in theii model, they concluded 
that a reduction of all rate measurement uncertainties to t|S9f 
would result in an overan uncertainty of ±.30^ for a given 
calculation of o2one reduction. They suggest, however, that a 
factor of two is a realistic estimate of uncertainty fur results 
Irom present models. However, this type of model is in its 
infancy, and it is dubious whether the data base for the 
reaction rate constants is adequate at present. In addition, the 
more fundamental question to be answered is whether a l-D 
model is adequate for an oione perturbation calculation. 

Tlie quality of the models is expected to improve with 
changes in the data base. In fact, several of the gaps which 
need to be filled were discussed in a recent review of 
stratospheric modeling with regard to aircraft (Ref. 31). Some 
of the more prominent items are simultaneous measurements 
of individual species to validate certain aspects of die l-D 
photochemical models, instruments to accurately measure 
water vapor in the stratosphere, and improvements in model- 
ing techniques to account for the true (nonlinear) nature of 
the atmosphere. 

C. Impact of the Stratospheric NOx Injection 

The impact of the NO^ injection on the stratosphere was 
calculat'd, in terms of -.lOj. by two methods. I'rom a 2-D 
model created by the Lawrence Livermore Laboratory 
(Ref. 33). the result for an S.ST lleet Hying at an altitude of 35 
lent was linearly scaled down to an injection magnitude of one 
SS abort flight per year. The second approach used was based 
on a formula from Johnston (Ref. 34). for .A(),/AN()^ 
following a 20-km injection; -,A()j/0, = (I ’5) ^SO^/SO^. 
The result from this latter approach is less reliable due to 
numerous changes in rate constant data since the formula was 
derived, and due to dilferetrces in injection altitude. 


« 

SS Injection; 

Injection' = 4.6 X 10^" molecules NOj 




a 


4.6 X 10“ 
3.2 X 10“ 


(-17.24%) 


=-2.48 X 10'*% 


Johnson Correlation (Ref. 27) 

”o7 ^ ?”no 

AO3 I ANO^ 

2 X 10'^ ^ 2 X lO** 

_ 1 (4.6X 10“)-:-(7.08X 10“)X(2X 10 '^) 

aO. - _ _ 

^ 2 X 10 ’ 


= -2.6 X 10* molecules 
= -1.3 X I ()■'•% 

The -AO3 calculated is ~5 X 10* (SST model) to —3 X 
10* molecules cm”^-s‘'. Tliis represents a change of —10'* 
to 10 •*%. This constitutes a change in IO3I in the strato- 
sphere about 4 orders of magnitude less than a change that 
would be detectable (~0.19f) (Ref. 6). Thus, for all practical 
purposes, it can be said that there are no impacts on 
stratospheric o/one from such a propellant dumping. Tliis is 
not unexpected when one considers that e.xisting SSTs emit an 
equal amount of NO, in less than 10 hours of cruise flight 
time. 

0. NOx Plume Chemistry in the Stratosphere 

The release of N2O4 into the stratosphere is equivalent to 
the release of twice the quantity of NOj as the will 

rapidly (in seconds) undergo thermal decomposition: 


LLL 2-D Model (Ref. 33) 

SST Fleet; 

Injeciion = 3.2 X lO''* molecules yr“' 
A()3 = - 17.247, 


NO, + NO, 

2 4 2 2 


' .Assuming one abort per year wiih tire N3O4 propeffant 
■f (ir one aborted mission. 
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Tlic NOj tlius produced will be photoly/ed (J ~~ 10' ^s*' 1 
within a lew minutcii to t'omi an equal concentration of nitric 
oxide and atomic oxygen. Hie nitric oxide reacts with o/one 
to regenerate NOj while the atomic oxygen recombines with 
the large concentration of molecular oxygen to form o/onc. 
Tlie overall result is that NO and NOj rapidly attain a 
transitory steady-state equilibrium, which is accompanied by a 
net production of o/one. 

N,o, -*: no, 

2 4 2 

NOj + hp - NO ♦ O 
O + Oj + M -* Oj + M 
NO + Oj - NOj + Oj 

After this mitial production of o/one, the standard NO^ 
catalyzed decomposition of odd oxygen occurs. Hiis decreases 
the odd oxygen content within the plume. However, this 
pe'turbation will be higlily liKali/ed and its impact on ihe 
total column abundance of o/one insignificant. Ciut/eii 
(Ref. 18) used a 2-1) model to demonstrate that the observed 
Oj perturbation that followed the injection of NO^, via a solar 
proton effect, in the upper stratosphere/lower mesosphere 
could be adcqiiaiely explained, thus lending some credence to 
our understanding of upper atmospheric o/one chemistry. 

Wliile the concentration of NO^ remains higli relative to the 
normal ambient levels, the conceniralions of numerous radi- 
cals deviate from their normal pliotochemically controlled 
equilibrium concentrations; e.g.. the OIK’rr) : HO, ratio will 
be controlled by the NO + HOj reaction rather than the 
O + OII. and O + HO, reactions. However, these short-term 
perturbations are of little significance as the NO^ plume 
rapidly disperses, and normal photochemical equilibrium is 
rapidly reestablished. 

E. Neutral Chemistry for a High-Altitude Injection 
of NO. 

Two of the three abort segments discussed in Section II 
occur in the altitude region of about 50 to 120 km, well above 
the stratosphere, encompassing the mesopheric and thenno- 
spheric regions of the atmosphere. As discussed in connection 
with plume dispersal, this material will be carried downward 
into the stratosphere at a rate that is relatively rapid, 
depending on the exact choice of eddy diffusion coefficients 
used for the calculation. Following t!m mixing process, the 
effects in the stratosphere on the orone layer will then hi’ 
identical in character, though not in magnitude, to those 
previously discussed for die case of a direct iniection into the 


stratosphere. There are. however, certain effects of the 
higli-altitude neutral chemistry that are of interest. One of 
these processes will have Ihe etfect of destroying significant 
fractions of the injected NO, so that much of it will rearch the 
stiatosplrere in an unreactive fonn. 

We first review the sequence of events following the release 
of N 2 O 4 in the high- -.Ititude 'egion. By far the most dominant 
mechanism for N 2 O 4 destruction would be collisional disstKia- 
tion. as previously concluded in the section on stratospheric 
rffects: 

N,(K + M-*2N(), +M (18) 

2 4 2 

Tlie rate constant for this reaction is 4 X 10"’ exp(-5Sb0/T) 
cm-*/s. so that for a total particle (M) density ol 10'* cm'* 
and a temperature of 220 K (which are typical for 80-km 
altitude), dissociation of the entire amount will (Kcur in ap- 
proximately four minutes. In all probability, this process will 
be accelerated somewhat by photodissociation, 

NjO^ +lu'-» 2 NOj (H>) 

with a 7-valuc that we estimate to be of the order of 10 " *s* ' . 

Similarly, the pholiHlissociation of NO, . 

NO, + hi- - NO + O (l.H 

which occurs with a 7-value of 10^^ s"'. will dissociate the 
NO, in about two minutes. Tlie reaction 

O + NOj-NO + 0, (12) 


will be even faster than the photodissociation, owing to the 
higli concentrations of atomic oxygen (10'” to iO' 'em ■') 
the higlier alti'udes, cotnbined with the very rapid rate 
constant for that reaction (‘J. I X 10*'- cm*/s). 

The preceding discussion shows that the N ,()4 wili lx- 
converted to NO on a time scale of minutes. Coincident with 
the process of plume dispersal, as discussed in Section II, the 
NO will pariicipate in several neutral chemistry reactions (in 
addition to those previously described), the main features ol 
which are; 

NO-Hu'-*N + 0 7-10'*s‘‘ (20; 
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N K)j -N0 + 0 * » 5.5 X 10*' 2 cxp(-3:20/T> 


CD 

N + N()-Nj+() A = H.2 X 10‘" exp<-4IO/T) (22) 

Tlicse processes are unique to the higli-altitude region because 
uf ( I) the rapid rate of NO photolysis and (2lthe presence of 
atomic nitrogen at significant concentrations (10^ to 10’’ 
cm'-*). The latter property of these regions is important 
because it provides a natural sink of NO at these altitudes. 
Tlius, the calculated lifetime of NO in the presence of 
10** cm'^ of N is about one day. Obviously this time is sliort 
compared to the time required for the NO to diffuse out of 
the region to lower altitudes, which is of the order of weeks 
(see Section II). It follows, therefore, that a high-altitude 
injection of NO will be much less effective, in terms of final 
NO^ concentrations achieved, than the corresponding injec- 
tions directly in the stratosphere. Siiwv w have already shown 
that the dinrt stratospheric iniections are not harmful, the 
input to the stratosphere of iniections at higher altitudes is 
even less effective. 

Another consequence of the presence of NO will be the net 
local production of odd atomic oxygen, as opposed to the 
catalytic destruction predominant at lower altitudes. The 
difference arises from the fact that at higli altitudes the Oj 
concentration is relatively low (so that Oj NO is slow), 
whereas both J and |N| are higli. Thus, the rate of production 
of atomic oxygen by the reactions written above exceeds ihe 
rate of catalytic destruction of otti' .).xygen by NOj^ . 

By way of summary, we have concluded that (I) N ,04 
injected at a higli altitude will be largely destroyed before it 
reaches the stratosphere, and (2) there will be a transitory 
localized net production of odd ox'-gen. Because of the 
temporary and localized nature of the production of atomic 
oxygen, no practical or harmful consequences are expected to 
occur. 

F. NOy Ion Chemistry 

The terminology used when discussing the upper atnur- 
sphere rellects the subject matter. IJiscussions of neutral 
chemistry refer to the 50 to ‘)0 km and ‘>0 to 120 km regions 
ol Ihe atmosphere as the mesosphere and thermosphere, 
rellecting the ambient temperature as the distinguishing 
characteristic. Overlapping this region is the ionosphere which 
starts at 60 km and extends upwaid. This retlects the elevated 
concentration of ions that first become significant at 60 km. 

Our purpose in this section is to summarize brietly the 
significant properties of the natural ionosphere, and to 
describe the perturbations to be expected as a result of the 


pro^v'llanl release. We will be concerned only with the 
D-region. which extends from 60 to ‘>0 km, and the fc-region. 
which extends from ^ to 1 20 km. 

In the D-region, ionization results almost exclusively from 
photoionization of naturally-iK'curring NO by the Lyman-a 
line at 1215 A. 

NO -H. 1 ^ -» NO* ♦ e (23) 

(7 = 5 X lO"’ at zeroopti<al depth) 

With an ionization potential of only 9.25 eV, corresponding to 
1340 A. NO is the only atmospheric constituent that cai; be 
ionized by this radiation. Sliorter wavelengths do not pene- 
trate to these levels. Although NO* is the major primary ion 
formed in the D-region, it is not the dominant ion at all 
altitudes in this region because other ions are pr.rduccd by- 
secondary, change transfer reactions leading mainly to 
hydrated protons. HdljO)*, and hydrates of NO*. The NO* 
tends to be more dominant near the top of the D-region, 
whereas the hydrated ions become relatively mure prevalent at 
lower altitudes (Ref. 35). 

Above 90 km, in the K-region. radiation of shorter wave- 
length is present, and phutoioni/at.on of Oj becomes the 
dominant priKess (however, phutoiunization of N, is also 
important): 

Oj +lu'(X< 1027 A)- Oj + e (24) 

Neverthelefs. NO* continues to be a major ion in the H-region, 
because it is formed rapidly by such charge exchanges as: 


Oj+NO-NO*+Oj (25) 

The steady-state concentration of any given ion, in eitl;e>' 
region, is determined by the balance between the rates ol 
production and the rates of destruction of that particular ion. 
Production mechanisms are direct photoionization by sunlight, 
and charge exchange from other ions. Tlie major loss media 
nisms are charge exchange to other species and (for poly- 
atomic ions) dissociative recombination with electrons: 


Oj+e-*0 + () (26) 
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Since (•hoioioni/alion of NO is (he major source of ion 
pririluction in (he l)-layer, (he ilensKy of charged par(icles will 
be increased in (ha( region during (he period ol'iime in which 
(he concen(ra(ion of (he dispersing NO cloud is comparable (o 
or greaier (han (he ambien( NO concen(ra(ion, which is in (he 
range of 10'^ (o 10" cm~^. Based on our plume dispersal 
calculadons, (hese conceni radons will be mainiained only for 
abou( one day following release. 

The increase in (he charged par(icle densiry can be 
calcula(ed from (he following s(eady-s(a(e equaiion: 

^ = y(p»io(oion»lNO| - a|NO*l (e) =0 

where |e| is (he elec(ron concen(ra(ion. y(pho(oion) is (he 
raie ci>efficien( for pho(oioni/.a(ion of NO, and a is (he ra(e 
coefficien( for dissocia(ive recombinarion of (he ion-elec(ron 
pair. Since (NO* | = jej . we have 



A( zero oprical dep(h,y(pho(oion) is 5 X 10 ’ s" ' • decreasing 
(o aboui lO*** s ' a( ftO km (Ref. 36). Tlie quamiiy a is 
8X lO-’cm^/sai 220 K (Ref. 35). 

In (he K-region, (he effec( of added NO in increasing (he 
ion dcnsiiy will be somewhai less impoi(an(. owing io (he fad 
(hu( pho(oioniza(ion of Oj and Nj. radier (han NO. is (he 
major souice of ion produciion in (ha( region. Relaiive 
pho(oioni/a(ion raies for NO compared io O, and Nj in (he 
t-region show iha( ambieni NO accounis for '7o or less of (he 
(o(al phoioioni/aiion raie (Ref. 35). Tlius. ( follows (hai a 
signil'icani increase in charge densiiy will iK'cur only when (he 
added NO is aboiK (wo orders of magniiude greaier (han (he 
amoieni NO, i.e., greaier (han 10’* (o I0'“cm'^. Tliis 
condilion does noi prevail for any signiilcani perirxl of i>me 
following release. 

NO in (he K-layer will also have (he effeci of depledng (he 
densides of o(her ions such as Oj and Nj, since, as already 
mendoned. these species undergo change-iransfer reaciions 
wiih NO: 

Nj+NO-Nj+NO* (27) 

To sumnari/.e. addiiion of NO to the ionospheric region 
will ( 1 ) resuK in increases in (he ion densiiy. and (2) change 
(he relaiive ionic composiiion when (he conceniraiion of 


added NO is greaier (han (hai of (he naiural adiiosphere. 
However, as slu'wn in Seciion II, (he la((er condiiion prevails 
only for a shori duK, and only on a localized basis. When fully 
mixed, (he added NO concern ral ions are inconsequenlially 
small compared lo (he ambieni levels. Al (he preseni lime, (he 
(ransiiory effecis occurring on a localized basis are no( known 
lo pose any (hreal lo almospheric slabilily, nor lo produce 
any effecis of a harmful or delelerious naiure. 

V. Ruorine in the Atmosphere 

A. Background 

Tliis chapler deals wiih (he almospheric injeclion of a 
propellani (ha( coniains fluorine, principally as F,. The 
magniiudes of (he periiibaiions lo (he siraiospheric ozone 
layer ind (he mesospheric/lhermospheric ion conienl are 
based on a derailed consideradon of lluorinc and oilier neulral 
and ion chemisiries. 

Over 1.6 X 10* metric tons (1.7 X 10* tons) of lluorinc 
were used in the manufacliirc of 6.9 X 10* metric tons (7.6 X 
10* (ons)ot chlorotluoromelhanes(CFMs) worldwide in 1973 
(Ref. 37). About 815? of this was manufactured for aerosol 
propellani uses that would bring about a short-term release of 
the CFMs to the atmosphere. Wliile chlorine from (his and 
other sources has beer extensively invcsiigalcd and modeled 
for its effect on (he chemistry of the siratosjihcrc. fluorine has 
been largely ignored. Preliminary models of stratospheric 
fluorine chemistry have been made by Rowland and Molina 
(Ref. 38) and Siolarski and Rundel (Ref. 3'^'rThvie are no 
detailed models that have attempted lo describe any pcrtinba 
lions to the ambient concentrations of the mesospl'.-re by the 
injection of halogen-containing species. 

Rowland and Molina point out .ha( F an l FO can react 
with odd oxygen |0(*P) and O3I in a manner analogous lo 
chlorine; 

F + 0, -FO + Oj (10)* 

FO-i^O-F + Oj (13) 

However, they einnhasizc the differences in rales of reaction 
and energetics of these two halogens, principally with regard 
to abstraction reactions of F with CH4 and Hj (H atom 
abstraction from H,0 was neglected.) On this basis they 
predict that, relative lo Cl, very little odd oxygen will hi 

*AII read, tin numbers in this chapter are consistent with Tables 9 and 
10. 
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destroyed before F is converted to HF. They furll'er point oat 
that tlie HF molecules formed are very strongly bonded (bond 
dissociation energy = 135 kcal mole'* ). and therefore are not 
subject to OH radical attack. Also. HF does not photo- 
dissiK'iate at a significant rate in the I MX) to 22JO A re»ion. 
Titus, they concluded that, although the stratospheric concen- 
trations of FX and CIX resulting from CFM dissociation will 
be comparable, the FO, reaction chain removing Oj can be 
disregarded relative to the CIO, chain. 

Stolarski and Rundel (Ref. 3Q) continued in this area by 
generating a simple but incomplete model of stratospheric 
fluorine chemistry based on the following reactions: 

F + Oj - FO + Oj 
F + Hj -» HF + H 
F + CH^ -HF + CHj 

HF + hr' - H + F 
FO + O - F + Oj 

FO + NO - F t NOj 
FO + Oj - F + 20. 

HF + 0('D)-*0H + F 

Using available experimental and estimated reaction rate data, 
they calculated that the catalytic efficiency of fluorine, for 
ozone destruction, is less than that for chlorine by at least four 
orders of magnitude. Hence, they concluded that catalytic 
ozone destructi n due to fluorine is negligible relative to that 
due to chlorine and further that, since hydrogen-bearing 
compounds are present at parts-per-million levels. HF forma- 
tion will not cause signillcant RH compound depletion at 
stratospheric FX concentrations below parts-per-billion levels. 

Recently, several groups of investigators have endeav.tred to 
measure present-day fluorine profiles in the stratosphere. 
Zander, et al.. (Ref. 40) measured hydrofluoric acid in the 
upper stratosphere, and reported a value for the mixing ratio 
of ~3 X 10"'° above 28 km. Based on their ground and 
balloon-borne observations, they estimate over 509^ of the 
atmospheric HF is above 25 km. Morz. et al.. (Ref. 41 ) utilized 


balloons and aircraft-borne filter samplers to determine fluor- 
ide levels during all four seasons at altitudes of 15 to 40 km. 
The average of the observed mixing ratio values (v/v) ranged 
from ~I.5X10*'' at 17 km, to ~l.l XIO"'® between 25 
and 37 km. The data showed considerable variability (a factor 
of 2) at all altitudes above 17 km. More recently. Farmer and 
Raper (Ref. 42) used the Jet Propulsion Laboratory Higli- 
Speed Interferometer in a balloon-borne gondola to measure 
HF and the ratio of HF to HCl. In the 14 to 38 km region, 
they obtained a ratio of 0.1 ±50%. The HF profile as a 
functii)ii of altitude is illustrated in Fig. 18 (Ref. 43). 



Fig. 18. Profile of the HF Mixing Ratio In the Atmosphere 


d. Approach 

Atmospheric modeling is immensely complex, involving 
such factors as horizontal and vertical transport, solar radia- 
tion. known and estimated constituent concentiations and 
fluxes, known and estimated reaction rates, and others. .Several 
such models have been generated for chlorine to predict the 
impacts of increased atmospheric CIX (CIX = HCl + ClONOj + 
CIO <• Cl) from CFMs and from the Space Shuttle exhaust 
effluent. In using such models, an amount of CIX is assumed 
to be injected, and the resultant change in ozone levels is 
predicted based on the efficiency with which the CIO, (Cl + 
CIO) catalytically destroys odd oxygen. 

To assess the effects of the fluorine injection, we have 
taken advantage of the analogy to chlorine. The first step was 
to tabulate the heats of formation of relevant species, a^ listed 
in Table 9. A comprehensive set of stratospheric reactions that 

obiginai.w^ 

OF POOR 
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TaM* 9. Haata of fonnation for kay apaciaa Tabla 10. flaactfona conakfarad for tha fHiorlna modal 


Species K) 


' 

19 

>2 

0 

HI 

-65.3 

lO 

26 

'«2 

3 

O 

59.6 

«2 

0 


34.1 

o‘d 

ID4.K 

II 

52.1 

"2 

U 

oil 

9.3 

IIOj 

5 

11,0 

-57.8 

II 20 J 

-32.6 

NO 

21.6 

NOj 

7.9 

NOj 

17 

NOI 

-15.4 

IINO 3 

-32.3 


33.2 

(■II 4 

-17.9 

ClIjO 

3.9 

HOI 

-23.5 

ClljOO 

■r6.7 


involve lluorine-containing species (FX = HF, F. FO, FOj. 
HOF. and NOF) was then assembled (Table 10). from which 
were chosen tlnrse reactions (Tabie 1 1 ) judged to be im- 
portant. Figure 1‘) summarises the resulting fluorine chemis- 
try. These reactions were, in turn, used to formulate two 
possible chemical models (1 and 11). whereby the injection of 
lluorine into the stratosphere would result in the destruction 
('f odd oxygen. A set of analytic expressions was then 
developed that describes the catalytic efficiency with which 
FO, destroys odd oxycen at steady state. Using the chlorine 
reactions listed in Table 1 2. an analogous analytic expression 
was then developed to describe the steady-state catalytic 


Kcacliun 

.5//. 

kcal 

mole"* 

Comment or l(H)inute 

•1 j ♦ hv - F ♦ F 

38 

Dominant initiation 
priKTcss 

! j *^ 0 -F + F 0 

-14.6 

a 

F j ♦ H - F ♦ HF 

-98.4 

a 

F j + NO - 1 NO + F 

-18.0 

PoMibly an important 
initiation proccu 

•F ♦Oj-FO + Oj 

-27.1 

Odd oxygen destruc- 
tion and 1 0 generation 

♦F + Oj + M - FOj + M 

-16 

.Major l Oj formation 
mechanism 

•F + H - III + II 

-32.2 

b 

•F + lljO-lll +OH 

-17.2 

b 

♦1 +CII 4 -HI +CH, 

-33.2 

b 

F ♦ lljOj - HF ♦ HOj 

-109.5 

a 

1 ♦ HOj -> III + Oj 

-89.3 

a 

1 + HNO 3 - HF + NO 3 

-35 

a 

HI ^ hi- - F + H 

136.4 

Unimportarit strativ 
spheric process: low 
photon flux at relevant 
wavelengths. Important 
alxive 80 km; Lyman-a 
photolyus 

•HI +0('l))- 1 +011 

- 11.2 

Dominant strato- 
spheric HI sink: regen- 
erates active lluorine 

HI +011-1 +II 2 O 

17.2 

Unimportant; slow reac- 
tion rate 

HI +011 tv " > 2)- F + HjO 

- 

Fxothermic; possibly an 
iniporlanl process 

10 + hi^ F + 0 

52.6 

probably unimportant 
relative to other FO 
mechanisms 

10 + 10 - 1 + F + Oj 

-14 

Unimportant at low FO, 
concentrations 

•lO + O-l +Oj 

- 66.6 

Loss mechunism foi odd 
oxygen 

FO + O 3 - 1 + 0, + Oj 

-14.1 

c 

10 + Oj - i Oj + Oj 

-57.1 

c 

•FO + NO - 1 + NOj 

-20.7 

d 

10 + H - HF + 0 

-83.8 

imporlani process in 
the upper mesosphere 

FO + Hj - HF + OH 

-90.0 

1 and « 

1 0 + HjO - HI + HOj 

-28.5 

d .ind e 

10 + ni4-HI +CII3O 

-69.5 

d and t 

FO + HjOj-HI +OH+Oj 

-49.4 

a 

FO + HOj - HI + O 3 

-62.2 

a 
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Tabt* 10 (eonM) 


stfBiMptMffIc fkiorin# mocM 


Reaction 

A//. 

kcal 

mole"’ 

Comment or footnote 

10 ♦ CH4 - HOF ♦ CHj 

♦ 1.6 

d and f 

FO+ HjO-»HOI +H 

♦2.6 

d and f 

FO ♦ H, -> HOF ♦ OH 

♦ 17.6 

d 

FO ♦ HjOj - HOF ♦ HOj 

-11.9 

d 

1 0 ♦ HOj - HOF ♦ Oj 

-54.5 

d 

F'Oj -♦- hv -• FO -*- 0 

82.6 

probably uriimpui- 
tant relative to other 
F'Oj loss mechanisms 

FOj -f 0 - FO ♦ Oj 

-36.6 

c 

FOj ♦ Oj - FO ♦ Oj Oj 

-II. 1 

c 

FOj ♦ Hj ^ + „Oj 

-63.3 

1 

FOj ♦ HjO - HF Oj + OH 

-1.2 

1 

FOj + CH 4 - HI ♦ Of jOO 

-43.7 

S 

FOj ♦ HjOj- HI + Oj ♦ HOj 

-30.7 

1 

FOj + HOj - HI + Oj + Oj 

-73.3 

( 

FOj ♦ HNOj - HI -• -Jj + NOj 

-19.0 

g 

FOj + Hj -> HF H Oj 

-15.9 

g 

FOj + CH 4 - HI CHj -• Oj 

-16.9 

g 

FOj -V no - FO + NOj 

9.3 

d and f 

FOj + NO - FNO + Oj 

-39.8 

( 

FOj + H - HI ♦ Oj 

-10.2 

Important process in 
the mesophere 


•Denotes an important reaction. 

®L'nimportant due to low nonfluorine reactant concentratio- 
*’Si)tnificant loss mechanism for active fluorine. 

‘-'Possibly sifinit'icant mechanism for odd oxygen destruction. 

•^Inhibits odd-oxygen destruction. 

'Possibly significant loss mechanism for active fluorine: (see Section V). 
Although A// appears to indicate that the reaction is endothermic and 
therefore slow, the magnitude of the unceitainties in the heats of for- 
mation of l-'O, l-Oj, and HOI- are such that the reactions could at least 
be thermoneutral thus possibly allowing the reactions to proceed at a 
significant rate. 

*Will decrease l-'O catalytic efficiency if FO2 loss mechanism is 
significant. 


efficiency for the CIO, desiruclion of odd oxygen. The 
analytic expressions for FO, and CIO, were then compared to 
estimate the relative efficiencies of the two odd-oxygen 
scavengers. Tlius. making use of explicit perturbation predic- 
tions from one-dimensional models for chlorine, the effect of 
fluorine can be assessed. 

Subsections C and Dshow the steady-state expressions that 
have been devised to describe the catalytic destruction of odd 
oxygen by chlorine and fluorine, respectively. These analytic 


Reaction Reaction No. 


Fj ♦ hv -• F ♦ F 

19) 

1 ♦ Oj - FO ♦ Oj 

110) 

1- ♦ Oj ♦ M - FOj ♦ M 

(11) 

F ♦ RH -» HI ♦ R 

(12) 

FO ♦ 0 - 1- ♦ Oj 

(13) 

FO ♦ NO - F ♦ NOj 

04) 

FO ♦ RH - III ♦ HO 

(15) 

FO ♦ Oj - F ♦ 20j 

(16) 

FO ♦ Oj- FOj ♦ Oj 

(17) 

FOj ♦ Oj - FO ♦ 20j 

(18) 

FOj ♦ 0 - FO ♦ Oj 

(19) 

FOj ♦ RH - HI ♦ ROj 

(20) 

HI ♦ 0('D) - 1- ♦ HO 
RH = HjO, CII4 . Hj 

(21) 



Fly. 19. Schematic diagram of the stratospheric fluorine model 
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Tabta U. RMCtton* and numbaring aystam uaad In 
atratoaphartc cMorina modal 


Reaction 

Reaction No. 

0 ♦ O 3 - no + Oj 

(I) 

CTO + 0 - a + Oj 

12 ) 

CTO + NO - a + NOj 

( 2 ) 

Cl + KH HCI + R 

(4) 

HO + OH - Cl + HjO 

15) 

no + NOj + M - CTONOj + .M 

( 6 ) 

ClONOj ♦ hv - CIO + NOj 

(7) 

Rtl = Hj. CTi.,. MOj 



expressions represent the catalytic efficiency after the systems 
have attained photochemical steady state. For the case of Fj. 
photiKhemical equilibrium is reached after the Fj has been 
completely photoly/ed (this takes ~16 hours after injection). 

Tlie plume cliemistry prior to the attainment of steady- 
state conditions is qualitatively described in Subsection G. We 
have not developed our own one-dimensional model, where 
photochemistry is coupled to parameterized transport in the 
vertical direction, because this type of model is inadequate to 
describe cither the short- or long-term perturbations caused by 
the injection of a pollutant in a plume. Therefore, we calculate 
tlie plume volume (F.\ concentration) as a function of time 
(see Section 11) and then estimate the steady state ozone 
perturbation at a scries of times after the initial injectitrn of 
the pollutant utilizing the analytical expressions. 

C. The Stratospheric Chlorine Model 

Atmospheric chlorine has natural as well as anthropogenic 
sources. Natural sources include sea spray, volcanoes, and such 
others as CHjCl. Anthropogenic sources include the Space 
Shuttle main engines and CFMs. Basic kinetic chlorine models 
incorporate all the reactions and photochemical processes that 
are considered important in the 0,, MO,. NO,, and CIO, 
systems. 

Althougli some 10 to 20 reactions involving chlorine- 
containing species are used in the iuodel calculations, (he basic 
reaction scheme can be described by the reactions listed in 
Table 12. Using this matrix of reactions, the catalytic- 
destruction of odd oxygen can be described by an analytic 
expression for the rate of odd-oxy^i-n destruction. -J(0 + 
O3 )/<//: 

-J{0 + 0.) 

- j- =(t,(Cl)( 03 ) + A.j(0)(CI0)- AjlNOMClO) 


(A) 


Steady-state expressions fur |CI|, |CIO|. (CIONOjI and 
|HCI| can be derived: 


*,(0H)(HC1) 

*/(RHT ' 

A-,(Cl)( 03 ) + y(C10N0j) 

" A ,(0) + Aj(NO) + A^(NOj ) ('^) 


ICIONO,!^^ 


A^(N0,)(C10‘ 

y, ■ 

A^ICIHRII) 

A 5 (OH) 


(III) 


(IV) 


Substitution of expressions I and II into equation (A) yields an 
analytic expression that describes the catalytic elTiciency of 
CIO, in terms of the concentration of the dominant chlorine- 
containing species. MCI. 

, 2A-, A- (OMHO.HHCl) A,(0) 

i=— X 

Jt * A.,(RH) A,(0) + A,(N0) 


(B) 


Equation (B) is a key result, and will be compared with 
analogous equations for lluorine to assess the fluorine effect 
on ozone. Tlie expression, (B), is the same whether or not 
reactions (6) and (7) are included in the mechanism. Tlie 
effect of including these two reactions is to reduce IHClj. 
However, the theoretical models show that the (CIONO, |/ 
|MCI| ratio is typically small. 

D. The Stratospheric Fluorine Models 

Wliereas considerable work has been done on stratospheric 
chlorine chemistry, and there is a consensus as to the principal 
aspects of the model, stratospheric fluorine chemistry' is much 
less certain. Reaction kinetics are known for very few 
reactions. In our approach, analytic expressions have been 
derived for two possible fluorine models that are consistent 
with the available experimental data. Tlicse are then evaluated 
individually. Two models were used to avoid the complexity 
that would result if all the possible reactions were incorpo- 
rated into a single model. 
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Fluorine may conceptually be partitioned into two states. 
HF and^ FO,(x = 0. I, 2). HF is very stable and unreactive, 
^u.cicas FO^ can catalytically destroy odd oxygen and hence 
ozone. 

Our models represent stratospheric fluorine neutral chemis- 
try using the most important reactions and the best reaction 
rate data available. The point of the models is to determine 
how much of the fluorine remains as FO, (destroying odd 
oxygens) once equilibrium is reached. 

Since the factor determining the magnitude of odd-oxygen 
destruction is the partitioning of F as MF and FO^, the key 
factors in generating the models are those reaction rates that 
affect either the hydrogen abstraction reactions (FO^ + RH -» 
HF + ROj) or those reactions that compete with hydrogen 
abstraction reactions. 

1. Fluorine model (I). A distinguishing feature of this, 
model is that the FOj radical is formed only via the F + O, + 
M reaction (II). and is not a product of the FO + Oj reaction. 
It is assumed that the dominant channel in the FO + Oj 
reaction results in the formation of atomic fluorine and 
molecular oxygen (16) rather than the FOj radical and 
molecular oxygen (17). Model II differs in this respect, as will 
be seen. Tlie present reaction system can be written: 


Fj + Im-* F + F 

(‘>) 

F + O 3 - FO + Oj 

( 10 ) 

F + 0, + M - FOj + M 

(ID 

F + RH-HF+ R 

( 12 ) 

0 + FO - F + Oj 

(13) 

NO + FO - NOj + F 

(14) 

RH + FO-HF + RO 

(l-S) 


O 3 ♦ FO -» F + 20j 

(16) 

FOj ♦ O 3 - FO + 20j 

(18) 

FOj + 0 - FO + Oj 

O')) 

FOj ♦ RH -* HF + ROj 

( 20 ) 

HF + 0( ' D) - F + OH 

( 21 ) 


It will be noted that there are three reactions (12, IS, and 
20) by which FO, is converted to the inactive reservoir, HF. 
Tlie only one of these three reactions definitely known to 
tK'cur is (12); reactions (IS) and (20) have been included for 
the sake of completeness because tliey are energetically 
feasible, althougli the rates are unknown (see Table 13). In our 
assessment, several possible values of the rate constant for 
reaction (IS) will be considered to illustrate the effect of this 
reaction in reducing the catalytic efficiency of fluorine. 
Reaction (20). which has a similar effect, will not be 
incorporated explicitly into the analytic expressions, since to 
do so would introduce needless complexity into the algebraic- 
relationships. This is equivalent to assuming tliat A'2u=0. 
which gives a worst case evaluation of the fluorine impact on 
ozone. That is, if is in fact finite, then fluorine will be 
even less efficient for ozone destruction. 

In our proposed reaction systems, the only mechanism by 
which the in.->i tive HF is converted to one of the active forms 
of FO, is reaction of HF with electronically excited oxygen 
atoms (O'D). Photolysis of HF in the stratosphere is an 
unimportant process, and thus is not considered in our model. 
However, photolysis becomes significant in the mesosphere 
and thermosphere (see later discussion). Tlie only additional 
process that miglit convert HF to atomic fluorine at a 
significant rate is: 

HI- + OH(v''>2)-*HjO-hF 

This reaction, and its implications on stratospheric fluorine 
chemistry, are discussed in Subsection F of tliis Section. 


Tlie net rate of odd oxygen removal for Mtxiel I can be written: 

-<f|Oj -1/10 + O3I 

- jT- = jr~ =^o<F)( 03 )"*. 3 <OXFO) + A,,(FO)( 03 ) 


A ,^(FO) (NO) + k , „(.POj ) (O 3 ) + A, ,(FOj ) (O) 
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Table 13. Expaflmantai valuaa for lluorino f act ion rala conatanta 


Reference 

Rate constant/cm^ molecule'* s'* 

Temperature, 

k 

Footnote 

H + h j - III ♦ II 




Albright. Doduvov, Lavrovtkaya, 
Morzov and Tal'roze. 1969 
(Ref. 44)» 

(2.0 t 0.16) X 10'*°exp(-1262 t lOO/T) 

294 to 565 


Rabideau, llecht and Lewii, 1972 
(Ref. 45) 

(4.2 1 0.3) X 10'* 2 

300 


Preferred value 

2.4 X I0'*“exp(-1262/T) 

194 to 565 

a 

1 j ♦ NO - NOI + 1- 




Rapp and Johnston. I960 (Ref. 46) 

1.0 X 10' *^ exp (-(757 t 50)/T) 

195 to 298 


Kim, Maclean, and Valance, 1972 
(Ref. 47) 

7.8 X 10'** 

-325 


Preferred value 

None 


h 

1 + O 3 - 10 ♦ Oj 




Wagner, Zetzsch and Warnatz, 
1972 (Ref. 48) 

2.8 X 10'** exp(-226/T) 

253 to 365 


r ♦ Oj + M - 1 Oj + M 




Zeizsch, 1973 (Ref. 49) 

5.2 X 10'*^ exp(656/T). M = He 

272 to 362 


Arutyunov, Popov and Chaikin. 

(7 + 2) X 10'** cm** moleculc'*s'* , M = He 

293 


1976 (Ref. 50) 

(1.4 1 4) X 10'** cm** molecule'* s'*, M = Nj 

293 


Chen, Trainor, Center and l•■yfc, 

(5.4 t 0.6) X 10'** cm* molecule'*s'* , M = He 

298 


1977 (Ref. 51) 

(1.5 * 0.3) X 10'** cm* molecule'*s'*, M = Oj 

298 



(2.5 t 0.5) X 10'*° cm* molccule'*s'*. M = HI- 

298 



(5.0 1 0.6) X 10'** cm* molecule'* s'*. M = Fj 

298 


Preferred value 

1.25 X 10'** cxp(656/T). M = Nj 

272 to .362 

c 

10 + 10 - 21 - + Oj 




Wagner, Zetzsch and Warnatz. 
1972 (Ref. 48) 

3.3 X 10'** 

298 


Clyne and Watson, 1974 (Ref. 52) 

(8.5 ± 2.8) X 10'** 

:'-8 


Preferred value 

1.5 X 10'*' 

298 

d 

1 + Hj - 111 + H 




llomann. Soloman. Warnatz. 

2.63 X 10'*°exp(-805/T) 

300 to 400 


Wagner and Zetzsch, 1970 
(Ref. 53)“ 

1.77 X 10'" 

298 


Dodonov, Lavrovskaya, .Morozov 
and Tal’roze. 1971 (Ref. 54)“ 

(3.0 ± 1.0) X 10'** 

293 


1 oon and RcM. 1971 (Ref. 55) 

8.05 X 10'** exp(-1243/T) 

253 to 348 


Bozzelli, 1972 (Ref. 56)“ 

(2.5 ± 0.7) X 10'** 

-297 


Kompa and Wanner. 1972 (Ref. 57) 

6.3 X 10'** 

298 
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TaM* IS (coMd) 


Reference 

Rate constant/cm^ molecule*' s*' 

Temperature, 

k 

Footnote 

Rabideau. Hechl and Lewis. 1972 
(Ref. 58)* 

(6.6 1 1.7) X 10*'^ 

293 


dyne. McKenny and Walker. 1973 
(Ref. 59) 

K + Hj - MK ♦ H 

(2.5 t 50T) X 10*" 

298 


Igoshin. Kulakov and Nikitin. 
1974 (Ref. 60) 

1.55 X IO"'®exp(-544/T) 
2.5 X I0“" 

195 to 296 
296 


Lam Thank My, Peron and 
Puget. 1974 (Ref. 61) 

2.3 X 10”'^ 

300 


Preferred value 

2.0 X 10"'° 

200 to 400 

e 

1 ♦ HjO - HI ♦ OH 




Zetzsch. 1971 (Ref. 62) 

2.0 X 10”‘‘ exp(-200/T) 

- 


F ♦ CH^ - HF ♦ CHj 




Fettis and Knox, 1964 (Ref. 63) 

2.0 X 10"'°exp(-6n9.fT) 



Foan and Reid. 1971 (Ref. 55) 

6.6 X 10"" exp(-930/T) 

253 to 3<‘8 


Wagner, Warnatz and Zetzsch, 1971 
(Ref. 64) 

5.5 X 10"'“exp(-579/T) 
7.9 X 10"" 

298 to 450 
298 


Kompa and Wanner, 1972 (Ref. 57) 

7.1 X 10"" 

298 


Clyne, McKenney and Walker, 1973 
(Ref. 59) 

6.0 X 10"" 

300 


Pollock and Jones, 1973 (Ref. 65) 

>1 X I0'“* 

298 


Preferred value 

7.5 X 10"" 

298 



®The two values reported for A:(300 K) are in fair agreement. Therefore, the preferred Arrhenius expression reflects the mean of the 
rate constants at 300 K. and the activation energy reported by Albright, et al. (Kef. 44). 


*'These two values differ by an order of magnitude at 300 K. and would probably be in even worse agreement at stratospheric tem- 
peratures. However, this reaction is only important in determining the rate of removal of 1- j (assuming that the rate of this reaction 
is comparable to 7,). not the catalytic efTiciency of KOu. 

••The preferred value is based upon three factors: (a) the average of the three determinations of A' (~300 K) 5.7 X 10"^^; (b) the 
ratio of the third body efficiencies reported by Arutyunov et al. (Ref. 50) Ajj j; Aj^^ = 2; and (<-) the assumption that the activation 
energy is not dependent upon the identity of the third body. fa(N2) “ i'alHe) = value reported by Zetzsch (Ref. 49). 

*1 Although the value of A(FO + TO) reported by Clyne and Watson (Ref. 52) was obtained in a more direct manner than that of 
Wagner et al. (Ref. 64) and as such is less susceptible to error due to the presence of complicating secondary reactions and thus 
would normally be preferred. The value to be recommended in this assessment is a weighted average of the two studies. 


'The value at 300 K seems to be fairly well established as (2.5 t 1.5) X 10"' ' ; however, the values reported for E/R are somewhat 
more scattered (544 to 1243 K). Therefore, it was assumed that 4 - 2 x I0~’® and "hen E/R was calculated to yield a value of 
2.5 X 10-" at 300 K. 

*^The preferred value was derived from the data of Wagner et al. <Rcf. 64). Kompa and Wanner (Ref. 57). Oyne et al. (Ref. 59). and 
Pollock and Jones (Ref. 65). 
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= *,o(F)( 03 ) + (FO) j*,3(0) + *,^( 03 )-*„(N 0 )} 
MF03){*,,(03) + *„(0){ 


(C) 


Using a steady-state analysis, which will be valid at long times (discussed in detail later), the following equalities can be derived; 


|F|. = 


*j,(HF)(0 'D) + (F0){A;,3(0) + A-,^(N0) + A.',^(0^ 
)t,„( 03 ) + *,,( 03 )(M) + A-,j(RH) 


(V) 


IF01..= 


it,„(F)(03) + A-,^(F02)(03) + A:,,(F03)(0) 
k , 3(0) + k , ^(NO) + ^^(Rll) + * , ^(03 ) ' 


(VI) 


it,,(Oj)(M)(F) 


•^°:1-%-,,(03) + A,,(0) 


(VII) 


IHF|„ = 


A,,(F)(RH) + A-,j(FO)(HF) 


(VIII) 


Ar,,( 0 'D) 


Substitution yields; 


-Jio + o,! a.,,( 0 'D)(HF)|a-,„( 03 ) + a-,,(0,)(M){ |:a-,3(oi + :a;,^( 03 ) + *,5(RM){ 

A,3(RH)|/t,„(03) + A-, ,(0,)(M)} +A, ,(RH)jA,3(0) + A,^(N0) + A, 5 (RH) + A, JO 3 )} 


(D) 


It is useful to reduce equation D to a simpler form which represents the worst case (maximum possible efficiency) for the 
destruction of odd oxygen by FO^. 

The Worst case; This occurs if the removal of FO by RH (RH = H-,0. CII4. H, , or any other II containing species) is slow. i.e.. 

A,,- 0 . 


-</(OJ 2A,,(0'n)(HF) jA-,„(03) + A, ,(0,)(M){ >1 

“ jF' " A^H) ^ + A■j^^N^^^•, J0,)[ * 

Important analogies between the chlorine and fluorine systems become apparent upon comparison of equations (B) and (F). 
Quantitative comparisons are made in Subsection F, to assess the absolute catalytic efficiency of fluorine for the destruction of 
ixld oxygen. 

2. Fluorine model (II). Tliis model differs from the previous one in that it is assumed that the dominant channel for the FO 
and 0, reaction is production of FO, + 0, (17) rather than F + 20, ( 16). Again, in the derivation of an analytic expression for 
-(/(O + 0^)Uir. we assume for the sake of simplicity that A,j|(FO, + Rll -» IIF + RO, )-“ 0. The reaction scheme can be written; 
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Fj + hv -♦ F ♦ F 


(9) 


F + O 3 -* FO ♦ 0 , 

( 10 ) 

F + 0, + M - FO, ♦ M 

(II) 

F + RH -* HF + R 

( 12 ) 

0 ♦ FO - F ♦ 0, 

(13) 

NO + FO - NO, + F 

(14) 

RH + FO -* HF + RO 

(15) 

O 3 + FO -» FO, + 0, 

(17) 

FO + 0 , -»F0 + 20, 
i i 2 

(18) 

FO, + 0 - FO + 0, 

(19) 

HF + 0( ' D) - OH + F 

( 21 ) 


The equations and steady state concentrations that describe this reaction scheme can be written: 


dt 


= A,„(F)(03) + (FO) j*,,(0) + A:,,(03)-A:,^(N0)} +(FOj) {*, ^( 03 ) + *, ,(0)} 


(F) 


[FI 


ss 


*3,(0' D) (HF) + (FO) I k , 3(0) + * , ^(N 0 )( 
^ 0 ^ 7 + *, , 7 ^) (M) + *,3(RH) 


(IX) 


_ * , „(F) (O 3 ) + * , ,(F03 ) (O 3 ) + * „(F03 ) (O) 
« *,3(0) + * , . (NO) + *, 5(RH) + *, ,(03 ) 


(X) 


|HFI„ 


*,,(F)(03i(M) ! *,,(03)(F0) 
*,«( 03 ) + *„( 0 ) 


*,3(F)(RH) + * 


,,(F 0 )(RH) 


(XI) 


*j,(0'D) 
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(XII) 


Therefore, by combining these equatioiu, the following expression for ~d{0^)/dt is obtained: 


-«f(0,) ^ W (HF) j A: , }2A:,3(0) f *,,(RH) t 2A:,,(03)| 
*,,(RH) j*,„(03) + A:,,(0,)(M)} +*,3(RH) |A:,3(0) + Ar,,(RH) + *,^(NO)| 

The worst case occurs when k^^ =^0, which simplifies to: 

-J(0^) 2ikj,(0'D)(HF)|A:,„(03) + k,, (Oj)(M){ ^ j*, 3(0) + *, ,(03)! 

~~dT " ^ |it,3(6) + A:,^(NO){ 


This equation is discussed in Subsection F in terms of the absolute catalytic efficiency of FO,. 


3. Summary of wont case equations for the fluorine models. This situation occurs when both A:, 5 (FO -t- RH) and *:o(FO + 
RH) are assumed to be zero. TIterefore, -J(0 + 0^)ldt can be written: 


Model (I): 


2A: j , (O' D) (HF) \k , ^( 03 ) + A , , (Oj ) (M )[ 
A, 3 (m 


|A,3(0) + A„(N0) + A„( 03 )| 


(E) 


Model (II): 


2 Aj,( 0 'D)(HF)}A,„( 03 )+A|,( 02 )(M)| 

”a, 3 (RH) 


|*7,(0)t*7;7Noif 


(I.) 


These expressions are similar, the only difference arising from the relative importance of reactions (16) and (17). Both of these 
reactions refer to the FO + O 3 reaction, but give different products (F and FOj), respectively. 

E. Evaluation of Rate Constanta for Key Fluorine 
and Chlorine Reactions 

Tlte rate constants that have been experimentally determined for the key fluorine and chlorine reactions are summarized in 
Tables 13 and 15 respectively. Table 14 summarizes our estimates for .'bose fluorine reaction rate constants that have not been 
determined experimentally. These rate constants are needed in Subsection F to evaluate the relative efficiencies of CIO, and FO,, 
which are predicted fro..; (he analytic expression derived in Subsections C and D. 

Notes accompany each table to explain our selection of the preferred values. 
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TaMt 14. Ertln m Ua of rat* oonatanto tor IhoM stratooptwrto 
fkiorino rMCttons lacking aiparknanlal data 


Rea 'tion 

Rate constant/cm^ 
molecule" *»■* 

Footnote 

0 FO - F -f Oj 

5 X 10"' ‘ 


• 

NO ♦ FO -« F' * NOj 

2x 10"" 


8 

O3 ♦ FO F ♦ 20j 

2x 10'" 

(i) 

b 


2x 10"'* 

(u) 



10"'^ 

(iii) 



10"'* 

(iv) 



0 

(v) 


O3 * FO - FOj Oj 

5 X 10'" 

(i) 

b 


2x lO'** 

(u) 



lO'*^ 

(iii) 



10-'* 

(iv) 



0 

(V) 


FO ♦ RH - HI- + RO 

0 

(i) 

c 


10'" 

(ii) 



10''* 

(III) 



lO'** 

(iv) 


FOj + O 3 - FO 20j 

5 X 10'" 

(i) 

d 


10-'2 

(ii) 



10-'* 

(iii) 



0 

(iv) 


FOj 0 - FO Oj 

5 X 10'" 


e 

FOj + RH - HI + ROj 

0 


f 

HF +0('l))-OH + F 

1 X 10-'° 


K 

HI + he - H 1 

J- lo-** 
080 km) 


h 

F'j -t^hi--! +1 

/ - 5 X 10'* 

(altitude independent) 

i 

H + FO - HI + 0 

5 X 10'" 


j 

H FOj - HI * Oj 

5 X 10'" 


k 


“These estimates arc probably accurate (2a) to witl.in a factor of 3 and 
are based upon the assumption that the reactivity of l O with 0(^P) 
and NO is similar to the reactivity of both CIO and BrO with 0(^P) 
and NO (provisional evperimental data supports this assumption). The 
analogous experimentally determined rate constants for CIO and UrO 
are: Aj(0 + CIO) = 5.0 X lO"". *3(NO + CIO) = 1.85 x lO”"; 
*(0 + BrO)=(3 1 lactor of 3) X lO"*'. and A:(NO + BrO) = 2.1 x 
10“*'. These rate constants are expressed in units of cm^ molc- 
-ulc~’s"’ and were measured at ~300 K. The temperature depen- 
dence of the rate constants for these radical-radical processes is 
expected to be small. 

*’The I'O + Oj reaction has two possible pathways which arc extv 
thermic, resulting the production of I- 20j (16) or I'Oj + Oj 


(17). Although this reaction has not been studied in a simple direct 
manner, two studies of complex chemical systems have in*'erred u>me 
kinetic information about it. Starrico et al., (Kef. 66) measured 
quantum yields for o/one destruction in I'j/f^j mixtures, and attrib- 
uted the high values, ~46(I0, to be due to the rapid regeneration of 
atomic fluorine via the l O + Oj -• I- ♦ 20j reaction (16). However, 
their results are probably also consistent with the chain propagation 
process being l O + l O -• 21' + Oj (the latter reaction has been studied 
twice by Wagner, et al., l972;Clync and Watson. (1974), but although 
the value of |l Iconsumed '» known to be close to 

unity, it has not been accurately determined). Consequently it is 
impossible to ascertain from the experimental results of Starrico et al.. 
whether or not the high quantum yields for ozone destruction should 
be attributed to the I'O + Oj reaction producing cither F ♦ 20^ or 
F'Oj * Oj (this prtscess is also a chain propagation step if the resulting 
I Oj radical preferentially reacts with ozone rather than with either 
FO or itself.). Wagner et al.. 1972. utilized a low-pressure discharge 
flow-mass spectrometric system to study the I + Oj and I'O ♦ FO 
reactions by directly monitoring the time history of the concentra- 
tions of F. FO. and Oj. They concluded that the FO + Oj reaction 
was unimportant in their system. However, their paper does not 
present enough information to warrant this conclusion. Indeed, their 
value of A:(F'0 + FO) of 3 X 10"* ’ is about a factor of four greater 
than that reported by dyne and Watson. 1974, which may possibly be 
attributed to either reactive impurities being present in their system, 
e.g.. 0(3/’). or that the FO + Oj reactions (16. 17) were not of negli- 
gible importance in their study. Consequently, it is not possible to 
determine a value of either A 16 or from existing experimental 
data, rherefore. in our assessment of the catalytic efticiency of FO^, 
a range of values for Aj^ and Aj.^ ssill be assumed (i through v). It is 
worth noting that the analogous CIO O3 reactions arc extremely 
slow (~l(r*** cm^ molecule" ’ s' * ). and only an upper limit of 9 x 
10"''* has been set for A(BrO ♦ O3). 

‘■There arc no experimental data for these reactions. Thus, it is difllcult 
to estimate the likely m.-’gnitude of these reaction rate constants. 
Some analogous CIO reactions have been studied and found to be very 
slow. e g.. CIO + Hj < 8 X 10"'* (670 K). CIO -t^CH^ < 4 X 10"'* 
(670 K). These values would correspond to upper limits of ~10'** 
and 10"''’ respectively at 230 K (stratospheric conditions). The 
fluorine reactions may be somewhat faster due to their greater exo- 
thermicity (due to the HF bond being stronger than the HCI bond), 
'riicrcl'ore in this assessment we have adopted a range of values for the 
10 RIl r.ite constants to test the sensitivity to these reactions. The 
case considered to be most important is that when A|j is assumed to 
be zero. 

‘^No experimental data. A reaction that has been studied that can be 
considered somewhat comparable is HOj -f O3. The literature values 
for this reaction are in rather poor agreement (especially the values of 
/■//?). but do agree that it is slow at 300 K with a value of -“1.5 X 
10"’*. A range of values can be used in the assessment. In the worst 
case (where Ajq = 0). it is irrelc ant what value of Aj ^ is used because 
(he only reaction considered (FOj Oy and I (), -*■ O) both yield the 
I O radical as a product. Only when A,q is given a I'inite value does the 
magnitude of A | ^ become important. 

®The rate constant for such a radical-atom process is expected to 
approach the gas collision frequency. As for A,g. the magnitude of 
A|^ is only important for nonzero values of Aiq. 

*^No experimental data for these or any analoeous reactions. 

^Assumed to be comparable to most other O’D reaction rate constants 
that approach the gas kinetic collision frequency. 




**Thi> J value hat been derived hated un a etiimated value fur ollll I at 
1 2 1.6 n*n leMrapolated from the ahturption tpectra uf MCI. MBr and 
HI), and the tlux of Lyman-a radiation. Thit it the dominant phutiv 
chemical lent mechanitm fur III in the upper meiotphere (>80 km). 

'Thit J value it bated on the abturption cruu-tectiunt reported fur I'j 
(Ref. 67). and the calculated J value of 2 X 10~'*i“* for CIj. Thit 7 
value it almoit independent uf altitude, at a major portion uf the 
abturption tpectrum it at waveleniitht greater than 3(X) nm. 


*Thete rale cunttanlt are expected to appiuach the fuit culliiiun fre- 
quency. Important luti procetiet fur active I'O^ in the metotphere. 


TaM* 15. Rala cofwUmts lor rooettona in atratoaphartc 
chlorina modal 

Reaction 

Rate conitant/cm^ molecule" 

O 2 + h*' -• Cl + Cl 

/■> 1 X lO"-* 

a ♦ O 3 - 00 + Oj 

2.7 X 10"" exp(-257n') 

CTO ♦ 0 - a ♦ Oj 

7.7 X 10"" exp(-I30/T) 

no ♦ NO - a + NOj 

1.0 X 10"" exp(-H200i 100)/T) 

0 ♦ CH 4 - Ha ♦ CH 3 

7.3 X 10"" exp(-1260/T) 

a ♦ HOj - HO + Oj 

3 X 10"" 

HC+OH-a + HjO 

3.0 X 10*" e\p(-»25/n 

no + NOj+ M - noNOj 

+ M. 3.3 X lO"*-* T"^ -’** 

M • Nj 

1 +8.7X 10 "’ * 

nONOj ♦ In; - no ♦ NO. 

2 y~ 3 X 10"* s"' (20km) 

~4.6 X 10"* (30 km) 

~ 1.9 X 10-*s"‘ (40 km) 
~4.2x 10"^ 5 *' (50 km) 

NOj + hv - NO ♦ 0 

7~ 5 X 10"* s*‘ 


F. Assessment of Fluorine Catalytic Efficiency 

Table 16 shows some typical values of the concentrations 
of atmospheric species that are predicted by the one- 
dimensional photochemical models, as a function of altitude. 
The values of these concentrations were combined with the 
preferred rate constants sliown in Tables 13 to 15 to 
numerically evaluate values of -</(0 + Oj)/Jf for both the 
fluorine (models I and ID and chlorine systems.' The calculated 
values of -c/(0+ 0_j)/Jf, which are shown in Table 17. 
represent the worst case situation where both A', j and are 
assumed to be zero. Tlie different values shown for -J{0 + 
Ojl/tff. (i) througli (v), correspond to a range of values for A, ^ 
(model I) and A , 7 (model ID as indicated in Table 14. Table 
18 and Fig. 20 show the ratio of the efficiencies of the 
fluorine and chlorine systems, p(F0^)/p(C10^)= 

( fluorine )/-t/( 0 ^)/i/f (chlorine) on a per molecule basis. 


From Tables 17 and IK ii can be seen that the values uf 
-<7(0, y<7/ calculated for model II are significantly greater than 
those fur model I. Tlie difference in catalytic efficiencies arises 
because in model I, F atoms are continually being recycled by 
reaction 16 (FO+ O 3 -* F ♦ 2 O 2 ). and are thus subject to 
being rapidly converted to HF via reaction 12. In model II, 
this is nut the case, because F atoms are nut involved in the 
cycle that destroys odd oxygen: 


FO 4 O 3 

- FOj + Oj 

117) 

FO 37 O 3 

-FO + 2 O 3 

(18) 

FO 3 +0 

- FO + Oj 

(I'D 


net: FOj + O (O 3 ) - FO + O, (20j) 

In both models (especially II) tl>e catalytic efficiency is 
greatest for high values of A(FO+ O 3 ). Table I9 shows how 
the catalytic efficiency of FO^ decreases (especially with 
model II) when the rate constant for tlie FU RH reaction 
(IS) is assumed to be nonzero. This reaction (15) provides an 
additional route by which active FO^ is converted to HF. A 
nonzero value for the rate constant of the FOj + RM reaction 
( 20 ) would have a similar effect to that of reictiun ( 15 ). 

Tables 17 and 18 show that even in the worst c se. the 
catalytic efficiency of CIO, normally exceeds that of FO,. 
Tliis difference can be attributed to the greater reactivity of 
atomic F. compared to atomic Cl. towards RH. and the inert 
nature of HF compared to HCI. Tlie only significant loss 
mechanism for HF is reaction with electronically excited O'D 
atoms, wliereas HCI is readily attacked by OH radicals, 
producing atomic chlorine. P i inert nature of HF. and the 
reactive nature of atomic F manifests itself in a larger 
percentage of FX being tied up in the form of HF compared to 
the percentage of CIX in the form of HCI. 

As it is impossible to estimate the most probable values for 
the rate constants of the FO + O 3 ( 1 6 , 17 ), FO + RH (15) and 
FOj + RH (20) reactions we shall assess the worst case, 
knowing that in all orobability the true catalytic efficiency of 
FO, (especially in n odel II) is at least a factor of 20 to 400 
lower. 

Table 20 summarizes the mixing ratios of FX that would 
result from a 50-s dump initiated at 43 km (this is referred to 
as case (a) in Section II). At times shorter than 10 days, the 
concentration of FX caused by propellant dumping is greater 
than ambient (see Fig. 3) where a mixing ratio of ~2.5 X 
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origin ALPAG^ 

OF POOR QUAU«» 


Spcciei 






'Vllitude. km 







20 

25 


30 



35 


40 

45 

M 

1.95 X I0“ 

8.4 X lO'’ 

3.76 X 10* * 

1.72 X 

lo'* 

8.2 X 10'* 

4.1 X 10'* 

Oj 

3.9 

X 10*’ 

1.7 X lO'* 

7.5 

X 

lO'* 

3.4 

X 

lo'* 

1.6 X lO'* 

8.2 X lO'* 

o'd 

5.6 

X 10*' 

3.0 

1.2 

X 

10 

3.8 

X 

10 

8.4 X 10 

1.2 X 10* 

O^f 

2.0 

X 10" 

1.4 X lo’ 

6.2 

X 

lO’ 

2.4 

X 

10* 

8.0 X 10* 

2.2 X lO’ 

HjO 

6.8 

X 10'* 

3.2 X 10'* 

1.5 

X 

10'* 

7.2 

«\ 

lo" 

3.9 X lO" 

2.1 X lO" 

CHa 

2.3 

X 10'* 

“.Ox 10 " 

3.0 

X 

10" 

1.1 

X 

to" 

4.5 X 10'® 

1.5 X lO'® 

»J 

1.3 

X 10'* 

6.3 X 10" 

2.8 

X 

10" 

1.2 

X 

lo" 

4.9 X lO'® 

2.0 X lO'® 

HO, 

4.9 

X 

o 

6.8 X 10* 

8.5 

X 

10* 

7.8 

X 

10* 

5.0 X 10* 

3.4 X 10* 

NO 

5.8 

X 10* 

7.3 X 10* 

9.8 

X 

10* 

9.7 

X 

10* 

9.5 X 10* 

6.7 X 10* 

o 

z 

1.9 

X 10® 

3.3 X lO’ 

3.6 

X 

10’ 

1.8 

X 

10’ 

4.8 X 10* 

6.0 X 10* 

OH 

1.1 

> 10* 

1.7 X lO" 

3.0 

• X 

10* 

6.3 

X 

10* 

1.4 X 10* 

1.6 X 10* 

Temp 

4.2 X 1C‘* 
218 

5.2 X lO'* 
223 

3.8 X 
229 

10'* 

1.8 X 
234 

lo'* 

5.4 X lo" 
248 

1.4 X lO" 
264 


10"'® (v/v) has been observed at 40 km. Recently (January 
1977), 2 NASA chlorofluorometi.ar.i assessment wmkshop 
was held to evaluate the stratospheric ozone perturbation 
cai.sed by the CFMs. Nine research groups calculated the 
steady-state ozone perturbation that would result from the 
continued release of CFMs at the 1 975 rate using one- 
dimensional photochemical models. The steady-state concen- 
trations of CIX produced from CFM photodissociation was —5 
ppb and the range of AOj that was predicted from the models 
was 5 to 97r. This assessment assumed that the fluorine 
released from the CFMs did not perturb the ozone concentra- 
tion. Since this workshop, there has been a redetermination of 
the NO + HOj -♦ NOj + OH rate constant that results in a 
doubling of the predicted catalytic efficiency of CIO,, i.e., 5 
ppb CIX reduces the integrated ozone column by 10 to 18%. 
For this assessment we shall assume that the percentage 
reduction in ozone is lineariy proportional to the CIX mixing 
ratio (this is not strictly true) and that 5 ppb of CIX causes an 
18% reduction in the integrated ozone column. Figure 21 
shows the results of the NASA CFM perturbation calculation, 
and illustrates how the ozone concentration is significantly 
reduced at 40 km, but is increased between 20 to 25 km 
predominantly due to the "self-healing" effect (because of the 
loss of ozone between 30 to 50 km, more ultraviolet radiation 
penetrates to tlie lower stratosphere, resulting in more 
low-altitude production of ozone). The fluorine situation is a 
little different from that of chlorine for several reasons: (a) at 
times less than a year, the fluorine resulting from a propellant 
dump is largely confined to the atmospheric layer in which it 
was initially deposited, whereas the chlorine resulting from 


CFM photodissociation is ali.iost uniformly distributed verti- 
cally (constant mixing ratio). For the sake ol simplicity, we 
have assumed that the fluorine is well-mixed vertically. This 
will tend to overestimate the artual ozone loss by a small 
factor ~(2 to 4); vb) It can be seen from Teble 17 that the 
magnitude of -t/(0,)/Jr in the chlorine model monotonically 
increases with altitude, whereas in fluorine model II the 
converse occurs. Consequently, it is difflcult to obtain an 
accurate value for the catalytic efficiency of FO, relative to 
CIO, when their relative eifciencies vary with altitude. In this 
as.sessment, the relative efficiency of FO, compared to CIO, is 
that predicted in Table 18 ai tiic mean altitude of the released 
cloud. Table 18 (model II) shows that p(FO,)/p(CIO, ) < 
2,1 (-3) at 40 km. After one day, when the plume diameter is 
~100 km, the ozone perturbation within the cloud will be 
p(FO,)/p(CIO,) X p(FO,)//j(CIO,) X 60% = 2.1 X 10'^ X 
2.l(-8)/5(-9) X 60% = 0.5%. The value of 60% fr.r ozone 
reduction at 40 km is taken from Fig. 21. However, the 
integrated ozone column will have been reduced by only 
«0.1 5% (due to the self-healing effect at lower altitudes, and 
because most of the ozone layer is between 20 and 30 km and 
the FO, cloud does not penetrate this region). Consequently, 
it can be seen that even the short-term localized perturbation 
is small. After one month, the FO, concentration due to the 
propellant dump is only ~5% of the background level, and the 
ozone perturbation is predicted (worst case) to be <10“*%. 
After five years the propellant is expected to be well-mixed, 
both horizontally and vertically, resulting in a mixing ratio of 
^10“*“* (more than 4 orders of magnitude less than present- 
day ambient). Propellant dumps initiated at sliglitly lower 
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TMt 17. CtlculHd valuM for tho chongo In odd oxygon 
gUgf ■ ilraliMptiwte Hunting ln|Mlton* 


AlIF 

(ude. Model 1 

km 

Model II 

CIO, model 

25 I** 2.4 X l(T* X |HF| 

4.2 X I 0 -*X (HFl 

9.2 X 10 ** X (HOI 

U 2.4 X lOr* 

1.6 X 10 ** 


iii 2.3 X 10 -* 

7.9 X 10 *’ 


Iv 7.2 X 10 *’ 

9 2 X 10 *’ 


v 1.2 X I 0 -* 

1.2 X lO*’ 


30 i 4.5 X I 0 -* 

3.6 X 10 “* 

8.4 X I 0 -* 

U 4.5 X I 0 -* 

1.5 X 10 ** 


Ui 4.3 X I 0 -* 

6.6 X 10 *’ 


iv 1.4 X 10 ‘* 

1.4 X 10 *® 


v 6.0 X 10 ’’ 

6.0 X 10 *’ 


35 i 7.1 X 10 -* 

2.0 X I 0 -* 

5.2 X 10-3 

u 7.1 X I 0 -* 

8.0 X 10 “® 


Ui 6.3 X 10 *® 

3.9 X lOr’ 


iv 2.9 X 10 -® 

2.9 X 10 -® 


V 2.7 X i(r® 

2.7 X lOr® 


40 i 6.7 X 10 * 

3.0 X 1 (T* 

1.4 X 10 ** 

U 6.7 X 10 *® 

1.3 X 10 ** 


ui 5.7 X 10 -® 

1.0 X ir'' 


iv 4.6 X 10 ‘® 

4.6 X nr® 


V 4.6 X ir® 

4.6 X 10 *® 


45 i 4.7 X 10 *® 

2.7 X 10 *® 

2.5 X 10 ** 

u 4.3 X I 0 -® 

1.5 X 10 *’ 


Ui 4.3 X I 0 -® 

4.7 X ir® 


iv 4.3 X 10 *® 

4.3 X 10 *® 


V 4.3 X 10 -® 

4.3 X lO*® 



^Calculated values for -d(0 ♦ Oj)ldt: worst cases (*| j • *20 " 
^These lower-case roman numerals correspond to the (i) - (v) range of 
values estimated for the FO Oj rate constant shown in Table 14. 

altitudes may cause slightly greater ozone reductions because 
p(FO, )/p(C10, ) is a factor of 200 greater at 30 km than at 40 
km. However, even if the catalytic efficiency of FO^ were 
equal to that of CIO,, the concentration of FO, resulting from 
a propellant dump would be so low that no significant ozone 
depletion would result after the plume had even partially 
dispersed. 

It is conceivable that, besides conversion of HF to active 
FO, via reaction with O' D, an additional process could be: 

HF + OH(v">2)-H20 + F 


iMIM IV* 1 119 fvPn PI COTWyilC XffNMflOMX Of Vw nUOfWlV 

■ — « — - . . 

ono emonno xywvfw 


Altitude, 

km 

Model 1 

p(FO,)/p<C10,) 

Model II 
p<FO,)/p(aO,) 

25 P 

2.6 X 10** 

4.6 

U 

2.6 X 10** 

0.17 

iii 

2.5 X 10** 

8.6 X 10*3 

iv 

7.8 X 10** 

1.0 X 10-* 

V 

1.3 X 10** 

1.3 X 10** 

30 i 

5.4 X 10** 

0.43 

U 

5.4 X 10** 

1.8 X 10 * 

iii 

5.1 X 10** 

7.9 X 10*® 

iv 

1.7 X 10** 

1.7 X 10** 

V 

7.1 X 10*® 

7.1 X 10*® 

35 i 

1.4 X 10** 

3.8 X 10** 

U 

1.4 X 10** 

1.5 X 10*3 

iii 

1.2 X 10** ' 

7.5 X 10** 

iv 

5.6 X 10*® 

5.6 X 10*® 

V 

5.2 X 10*® 

5.2 X 10*® 

40 i 

4.8 X 10*® 

2.1 X 10'* 

ii 

4.8 X 10*® 

9.3 X 10** 

lit 

4.1 X 10*® 

7.1 X 10*® 

iv 

3.3 X ir® 

3.3 X 10*® 

V 

3.3 X 10*® 

3.3 X 10*® 

45 i 

1.9 X 10*® 

1.1 X 10~® 

U 

1.7 X 10*® 

6.0 X 10*® 

iii 

1.7 X 10*® 

1.9 X 10*® 

iv 

1.7 X 10*® 

1.7 X 10*® 

V 

1.7 X 10*® 

1.7 X 10*® 


•These lower-case roman numerals correspond to ihe(i)- 
(v) range of values estimated for the FO -t- O 3 rate 
constant shown in Table 14. 

However, the maximum likely rate constant for the above 
process is and at that rate, combined with the 

predicted values for |OH, v" > 2| reported by Nagy et al. 
(Ref. 68) of ~10^ cm'^, the overall catalytic efiiciency of 
FO, would only increase by less than a factor of ten. 
Consequently, the catalytic efficiency of FO, would still be 
too low to be of importance in terms of releasing small 
quantities of F 2 from the Space Shuttle. 
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lO"* 10** 10*^ 10** 10** 10*' 10® lo' 


f(fO ' XCIO ) 

A M 

Rg. 20. Ratio of tha catalytic afficlenciaa of atratoaplMric 
fluorina and chlorina. 

G. Plume Chemittry 

Although the equilibrium eventually reached will resemble 
situations previously discussed, the propellant ejected from the 
Space Shuttle will initially exhibit "plume” chemistry, involv- 
ing localized and severe short-term perturbations of the 
ambient concentrations of the trace constituents. In tlie 
immediate vicinity of the Space Shuttle, a transient condition 
of a vcry-higli concentration of Fj will exist. Tlie time taken 
to reach steady-state conditions is governed by the rate of loss 
of molecular fluorine. Tlie dominant loss mechanism for 
molecular fluorine is photolysis, with a J value of ~5 X 
I0~®s"', with reaction with NO being a second minor loss 
process: 

Fj + Im -► F + F 


Fj + NO - F + FNO 
FNO + hi^ - F + NO 

Consequently, the time taken for 95% of the molecular 
fluorine ‘o be dissociated is about 16 h, assuming that 
photolysis is the only removal mechinism. It is dilTicult to 
ascertain the importance of the NO + fj reaction, as the rate 
constant is uncertain by at least one OiJer of magnitude (Table 
13), and, in addition, the NO in the plume may be 
significantly depleted if FNO photolysis is slow. However, if 
the NO + F 2 is rapid enougli to be competitive witli 
photolysis, its only importance is that the molecular flu urine is 
removed more rapidly and FO^ attains steady-state conditions 
earlier. Atomic fluorine formed from the processes shown 
above reacts with oxygen-containing species (Oj, O 2 . and O) 
to destroy odd oxygen, and hydrogen-containing compounds 
(RH = H 2 O, Cil 4 , H 2 ) to produce HF and localized high 
concentrations of odd-hydrogen radicals such as OH, CHj, and 
H. Tliese localized concentrations of odd-hydrogen radicals 
may far exceed the ambient levels (by many orders of 
magnitude ), and the plume will exhibit a chemistry dominated 
by radical-radical recombi.iation processes. 

While the plume is still higlily localized (Hiall a d:iy) the 
concentration of RlKHjO, C'H^, and H 2 ) will be almost totally 
depleted by the F + RH reaction, resulting in an effective 
increased catalytic efficiency for FOj,. Tlie ozone level in the 
plume will also be significantly, if not totally, depleted 
However, it is not worth carefully modelling this localizeu 
:-ituation as the condition is transitory, and the atmosphere 
rapidly recovers as the plume disperses. The injection of 
fluorine is in a region of the stratosphere where the chemical 
composition is controlled by photochemistry; consequently. 


Tf* 19. Changes in FO. catalytic atTiclancy as a function of 
nunzsfo rats constants for rsaction (15) (FO RH — HF -t- RO) 


Altitude, 

km 

*‘is 

*‘l6,17 - 
Model 1 

5 X 10"" 
Model il 

'‘l6.n ^ 
Model 1 

2 X 10"'* 
Model 11 

30 

0 

5.4 X 10"* 

0.43 

5.4 X 10"* 

1.8 X 10"* 


10-'* 

5.4 X 10"* 

2.1 X 10"* 

5.1 X ur* 

8.1 X 10"* 


10"" 

4.5 X 10*5 

2.3 X U)-'* 

1.8 X 10"* 

2.1 X 10"* 

40 

0 

4.8 X 10"* 

2.1 X 10"* 

4.8 X 10"® 

9.3 X 10"* 


10"'* 

4.8 X 10"* 

8.6 X lO""* 

4.4 X 10"* 

3.8 X 10"* 


10"" 

3.8 X 10"‘ 

1.6 X 10"* 

1.7 X 10"* 

l.V V 10"* 
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TaM* 20. Miiing ratioa of FX following a SO-a duration a)actk>n 
bafwaan 43 km and 30 km* 


Time 

|IX| 

h'X (uM 

% I'j (diiiociatcd) 

1 h 

5.0 X 10'* 

6.3 X 10~* 

17 

5.6 h 

7.1 X 10'" 

8.9 X 10"'' 

64 

1 day 

1.7 X 10’ 

2.1 X 10"* 

99 

10 day 

9.8 X 10* 

1.2 X I0*‘° 

100 

1 mo 

1.2 X 10* 

1.4 X 10“" 

100 

1 yr 

4.7 X lO'* 

5.9 X 10“*^ 

100 


*C'asc (a) SO-5 dump initialed at 43 km. 



AO/O3. % 


Fk). 21. Perturbation in the ozone concentration resulting from 
CFMb in the atmosphere. 

tliere will be a rapid regeneration of photodiemically pro- 
duced species such as Oj. 

H. Neutral Fluorine Chemistry in the Mesosphere 

Neutral fluorine chemistry in the mesosphere is similar to 
that previously described for the stratosphere. The major 
differences are that HF will undergo photolytic decomposition 
by short wavelength UV liglit (Lymcn-a at 1 21 .6 nm) at a rate 
much greater than its reaction with electronically e.xcited 
oxygen atoms (O'D). and ihat the concentrations of atomic 
hydrogen, nitrogen, and oxygen are much greater in tlio 
mesosphere than in the stratosphere. Consequently, this 
affects the relative importance of sjine of the reactions. The 
chemical scheme can be written: 


Fj + hi' *• F + F 

( 9 ) 

Fj + li - HF + F 

(22) 

F + O3 - FO + Oj 

(10) 

+ Oj + M - FOj + M 

(11) 

F + RH-*HF + R 

(12) 

0 + FO -* F + Oj 

(1.^) 

NO + FO - NOj + F 

( 14 ) 

Rll + FO - HF + RO 

( 15 ) 

FO + O3 -► F + Oj + Oj 

( 16 ) 

FO + O3 -* FO3 + Oj 

( 17 ) 

ll + FO-»HF + 0 

( 24 ) 

FOj + O3 - FO + 20 j 

( 18 ) 

FO3 + 0 -» FO + O3 

( 19 ) 

FO3 + RH -»HF+ RO3 

(20) 

hf + o'd-oh + f 

(21) 

llF + hi' -*H + F 

( 23 ) 


The Table 21 shows some typical values expected for 
concentrations of neutral constituents in the mesosphere (Ref. 
69). 

A careful analysis of the fluorine reaction scheme listed above 
shows that about Five of these reactions dominate and are 
sufficient to describe neutral fluorine chemistry. Pltotolysis 
remains the major loss process for Fj until UX) km, when 
re. ction with atomic H becomes comparable. Tlierefore, the 
half life of F^ in the mesosphere is comparable to that in the 
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k,3(0)(F0) + J„(HF) 
k, (03) + k,, (RH) 


Tabto21. ConoanlratloiM ol iwulral oonstHuania in 
motMUtM/Cffl* 


Altitude, km 

Spectei 

60 

80 

100 

|M| 

7.2 X to'* 

4.2 X 10'* 

1.1 X 10 '* 

|Oj| 

1.4 X 10* * 

1.1 X 10'* 

2.2 X 10*2 

( 0 'D| 

~ 3 X 10^ 

~1 X 10* 

1.2 X 10* 

103 p| 

7 X 10 ’ 

6X 10® 

4 X 10 " 

[O3I 

3 X 10’ 

1 X 10* 

1 X 10 ’ 

ICH4I 

3.2 X 10* 

6.3 X 10* 

- 

IH2I 

3.6 X 10® 

1.7 X 10® 

1.1 X 10 ’ 

|H.O| 

3.5 X 10*0 

5.9 X 10* 

7.7 X 10* 

(H) 

1.4 X 10* 

3 X 10 ’ 

1.1 X 10* 

(OH) 

3 X 10* 

1 X 10* 

-- 

(HO2I 

3 X 10* 

1 X 10* 

- 

(NO) 

5 X 10’ 

5 X lO’* 

5 X 10’ 

T,K 

253 

177 

210 


(F| 


SI 


IFO).. 


k,3(0) 


IHFl.. 


k,, (F)(RH) 


Substituting values for the rate constants and concentrations 
in the above equations shows that HF:F:FO at 80 km is 
1 :10"‘’:I0"*; i.e., HF is the dominant fluorine species, as was 
the case in the stratosphere. The quantity of Fj ejected into 
the mesosphere is too low to significantly perturb the 
concentrations of any of the neutral species after a day or so. 
For a few hours, the concentrations of RH and O 3 may be 
depleted in the plume, but, as the plume disperses, the 
atmosphere rapidly recovers. Consequently, we conclude that 
the only fate of F 2 injected into the mesosphere is downward 
transport into the stratosphere in the form of HF, with no 
adverse effect on the neutral contituents of the mesosphere. 


stratosphere, h. Due to the decrease in total gas density 
with altitude, the F + O 3 reaction dominates tne F + Oj + A/ 
reaction so that primary formation of FOj is not particularly 
important in the mesosphere. The fate of the FO radicals 
formed in reaction ( 10 ) will be dominated by reaction with 
O^P. so that reactions (14), (15), (16), and (17) can be 
neglected. Reaction (24) is of little significance due to the high 
[O^Pl/(H) ratio. Reactions (18) and (19) can be neglected 
because the rate of formation of FOj via reactions (II) and 
(17) is unimportant. Therefore, the reaction scheme can be 
reduced to; 


Fj + hn -» F + F 

(9) 

F + O 3 -FO + Oj 

( 10 ) 

F -I- RH - HF + R 

( 11 ) 

FO + 0 - F + Oj 

(13) 


I. Fluorine Chemistry in the Ionosphere 

In Section 11, it was sliown that release of N 2 O 4 propellant 
in the ionosphere would result in temporary, localized changes 
in charge density and ionic composition, owing mainly to the 
fact that NO (which is the major product) has a very low 
ionization potential. For F 2 release, however, no such effects 
are expected to occur because of the much higher ionization 
potentials involved. Tlie relevant fluorine species, along with 
the corresponding ionization potentials, are; 


Species 

Ionization potential. eV 


15.7 

F 

17.4 

HF 

15.8 

FO 

13 to 14 (uncertain) 

FOj 

12.6 


HF + Im -► F + H (22) 

The following equations can be derived for [F],,, {FOjjj, 
and (HF1„. 


Thus, there will be no photoionization of the fluorine species 
in the D-layer, because radiation of sufficiently short wave- 
length is not present. In the E ^gion, these species will 
probably be photoionized withy-values similar to those of N 2 , 
since the ionization potentials are similar. However, since the 
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fluorine concentrations will be small, the resulting effects on 
the charge density or ion composition will be quite negligible. 
Furthermore, the principal form in which the fluorine v/ill be 
present, HF, cannot be ionized by charge exchange from any 
of the abundant ions in the D- or E-regions (Oj^, Nj^, or 
NO* ), because the ionization energy of HF is greater than that 
of any of those ions. 

We conclude, therefore, that release of Fj in the ionosphere 
will not significantly perturb that region. 


VI. Summary and Conclusions 

If a Space Shuttle flight must be aborted before attaining 
escape velocity, the propellant for the payload would be 
ejected into the stratosphere or the ionosphere (which includes 
the mesosphere and the thermosphere). Tlie payload propel- 
lants may contain up to '1122 kg (6,0(X) lb) of Fj or up to 
3402 kg (7.500 lb) of N 2 O 4 . The present report is an 
evaluation of the effects of these injections on the atmosphere. 

The addition of NO and NOj to the stratosphere has 
previously been extensively investigated in the course of 
evaluating the effects of SST exhaust on the atmosphere. In 
particular, there is concern about stratospheric ozone levels 
and the assiK'iated shielding of harmful U\ -B radiation. 
Assuming a linear relationship between the stratospheric load- 
ing of NOj and the magnitude of the ozone perturbation, we 


have calculated the change in ozone expected to result from 
the space shuttle ejection of NjO^, based on the ozone change 
that is predicted for the (much greater) NO, input that would 
accompany large-scale operations of SSls. These calculations 
show that the effect on ozone is negligibly small. 

The N 2 O 4 may also be released in the ionosphere. The 
resulting effects on ion density and composition are discussed. 
Because of the localized and transient nature of theeffccts.it 
is concluded that this will result it. no adverse environmental 
impacts. 

Although no previous extensive efforts have been made to 
evaluate fluorine additions to the atmosphere, other halogens 
have been considered. In particular, large-scale efforts have 
gone into the question of chlorine-catalyzed reductions in 
stratospheric ozone. We have critically reviewed possible 
stratospheric fluorine reactions to evaluate the magnitude of 
fluorine-induced ozone destruction relative to the reduction 
that would be caused by addition of an equal amount of 
chlorine. Thus, using the magnitude of ozone decrease 
predicted for chlorine, the decrease in ozone that would result 
from the fluorine addition has been calculated. Because of the 
chemical properties of fluorine, and the low concentrations 
involved, the predicted effect on stratospheric ozone is 
vanishingly small. 

A similar evaluation was made for an ionosphere injection. 
No adverse environmental impacts are predicted. 
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